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ABSTRACT 


Three  metrics  for  predicting  the  time  required  to  identify  targets  in 
high -resolution  radar  pictures  were  developed.  One  metric,  based 
on  four  automatically  measured  variables  related  to  transmissivity 
of  positive  transparencies,  was  tested.  Through  multiple -regres sion 
analysis,  a  correlation  of  0.69  between  observed  and  predicted  iden¬ 
tification  times  was  obtained.  When  these  relationships  were  applied 
to  a  new  set  of  radar  pictures  and  new  test  subjects,  the  correlation 
coefficient  was  too  low  to  be  significantly  different  from  zero  with  the 
number  of  pictures  used.  A  principal  reason  for  the  low  correlation 
was  the  unexpected  correlation  between  the  four  transmissivity  vari¬ 
ables.  The  metric  also  had  an  insignificant  correlation  coefficient 
when  applied  to  optical  and  infrared  photographs  (0.07  and  0.04,  re¬ 
spectively).  It  is  hypothesized  that,  if  additional  variables  are  utilized, 
the  metric  examined  in  this  study  may  be  useful  as  a  base  to  develop  a 
more  effective  prediction  equation. 
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SECTION  I 


INTRODUCTION 


The  relationship  of  certain  physical  characteristics  of  radar  pictures  to  the 
time  required  for  humans  to  identify  designated  targets  in  these  pictures  was 
investigated.  As  used  here,  "target"  means  the  radar  return  from  a  single 
object  or  combination  of  objects  making  up  a  small  part  of  a  relatively  large 
radar  picture  containing  many  returns.  The  purpose  was  to  express  the  re¬ 
lationship  in  a  metric  that  would  maJce  it  possible  to  predict  target-identifica¬ 
tion  times  for  similar  radar  pictures.  A  secondary  purpose  was  to  determine 
if  the  metric  could  be  used  to  predict  the  time  required  for  identification  of 
targets  in  infrared  and  optical  photographs. 

This  was  an  exploratory  effort  supporting  the  continuing  research  into  the 
larger  field  of  image  interpretation.  Eventually,  establishment  of  the  com¬ 
plete  psychophysical  relationship  between  quantitative  measurements  of  pic¬ 
ture  content  and  human  response  could  aid  in  the  planning  of  reconnais sance - 
strike  missions,  and  could  contribute  to  the  development  of  automatic  equipment 
for  identifying  targets.  However,  the  exploratory  effort  reported  here  was 
limited  essentially  to  consideration  of  only  the  relative  lightness  and  darkness, 
detail,  and  size  of  targets  and  their  backgrounds.  It  did  not  include  such  target 
features  as  shapes  and  patterns  or  inferential  cues  that  would  bn  present  in  an 
operational  reconnaissance-strike  situation. 

The  departure  from  a  more  realistic  situation  was  influenced  in  part  by  a  re¬ 
quirement  that  the  selected  characteristics  be  measured  automatically.  To  a 
large  extent,  these  characteristics  were  chosen  to  reflect  the  stimulus  com¬ 
plexity  of  the  radar  pictures.  In  this  way,  the  experiment  was  not  confounded 
by  parameters  of  lesser  interest  to  this  particular  study,  or  by  those  that  could 
not  be  measured  with  state-of-the-art  equipment. 

The  emphasis  was  on  high-resolution  ( side -looking)  radar  pictures.  All  had  the 
same  general  resolution,  graininess,  acutance,  etc.  The  experiment  consisted 
of  the  following  general  steps; 

Measurement  of  human  target-identification  times  for  a  series 
of  high-resolution  radar  pictures  (section  II) 

Measurement  of  certain  physical  characteristics  of  these  pic¬ 
tures  (section  UI) 

Relating  these  two  sets  of  data  mathematically  in  a  metric  to 
permit  prediction  of  tar  get  -  recognition  times  for  similar  pic¬ 
tures  (section  IV) 

Checking  the  validity  of  this  metric  with  new  pictures  and  new 
subjects  (section  V) 

Checking  the  applicability  of  the  metric  to  infrared  and  optical 
photographs  (section  V) 

The  task  required  construction  of  a  combined  slide  projector  and  viewer  and 
of  a  flying  spot  scanner,  programming  of  an  analog  computer  tc  measure  the 
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variations  of  tiie  selected  physical  characteristics  of  the  pictures,  and  use  of 
a  digital  computer  for  data  analysis. 

Three  prediction  metrics  were  developed;  only  one  was  evaluated. 
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SECTION  U 


MEASUREMENT  OF  HUMAN  PERFORMANCE 


INTRODUCTION 

The  human  response  data  required  for  generation  of  the  prediction  metric  were 
recorded  in  identical  experiments  with  30  Strategic  Air  Command  bomb-navi- 
gators,  as  discussed  below.  The  other  data  required  for  the  metric  -  measure¬ 
ments  of  physical  characteristics  of  the  pictures  -  are  discussed  in  section  III. 
To  obtain  sufficient  subjects  to  make  the  data  significant,  the  human  response 
experiments  were  conducted  at  two  SAC  bases,  Griffiss  in  New  York  (11  sub¬ 
jects)  and  Pease  in  New  Hampshire  (19  subjects). 

The  experiments  were  somewhat  artificial,  in  that  they  were  neither  recon¬ 
naissance-strike  nor  image-interpretation  tusks.  In  an  operational  situation, 
more  information  about  a  target,  its  surroundings,  and  its  location  would  be 
given  to  the  subject  and  thus  he  could  use  associational  and  inferential  cues. 

For  example,  if  a  subject  knew  that  the  target  were  a  dam,  he  would  infer  the 
presence  of  a  lake  and  roads.  However,  the  experiments  did  involve  some  of 
the  perceptual  factors  inherent  in  reconnais  sance  -  strike  and  image -interpreta - 
tion  tasks,  such  as  target  size  and  detail,  and  contrast  with  the  background. 

GENERAL  PROCEDURE 

The  slide  projector -viewer  built  under  the  contract  was  transported  first  to 
Griffiss  and  then  to  Pease.  The  subjects  were  trained  to  bring  all  to  the  same 
level  of  familiarity  with  high -resolution  radar  pictures  to  prevent  contamination 
of  the  scores  by  learning  acquired  during  the  experiments.  Each  subject  was 
trained  individually  for  1  hour.  High-resolution  and  ppi  radar  pictures  of  the 
same  areas  recorded  simultaneously  by  the  same  airplane  on  Air  Force  Quick- 
Check  flights  were  used.  Operational  navigation  charts  of  the  same  areas 
aided  the  interpretations.  The  increased  resolution  and  different  reflectance 
properties  of  the  high- resolution  radar  pictures  were  emphasized.  The  dif¬ 
ferences  between  high-resolution  and  ppi  displays  of  the  following  were  cited: 
large  cities,  land-water  boundaries,  communication  and  transportation  links 
in  various  terrains,  geologic  features  and  drainage  patterns  in  mountainous 
areas,  and  small  towns  in  flat  and  hilly  terrain  and  in  terrain  confounded  by 
the  texture  of  vegetation.  High-resolution  radar  pictures  similar  to  those  to 
be  used  in  the  tests  were  shown,  and  typical  targets  were  pointed  out. 

Thirty  targets  and  the  30  high  -  resolution  pictures  from  which  they  had  been 
extracted  were  used  in  the  tests.  These  pictures  also  were  obtained  in  Quick- 
Check  flights.  They  were  taken  from  20,  000-ft  altitude;  each  covered  an  18- 
by  1  8 -nautical -mile  ground  area.  Seven  other  similar  targets  and  pictures 
were  supplied  for  practice  trials.  All  the  test  targets  and  pictures  containing 
them  were  prepared  on  positive  slides  for  use  with  the  projector -viewer . 

EXTRACTION  OF  TARGETS 

One  entity  in  each  picture  was  designated  a  target.  These  were  limited  to  air¬ 
fields,  industrial  complexes,  oil  tank  "farms,  "  small  towns,  and  military 
camps.  The  targets  had  the  same  size,  pattern,  orientation,  and  approximate 
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density  as  in  the  complete  pictures.  As  such,  they  were  pictorial  briefing 
cues  requiring  no  mental  adjustments  for  locating  them  in  the  pictures. 

Human -dynam ic  and  cybernetic  studies  suggest  that  a  human  is  at  his  best 
when  he  functions  in  this  way  as  a  simple  comparator  (reference  1). 

The  targets  were  prepared  by  projecting  a  complete  scene  and  tracing  the 
selected  target  image  on  white  paper.  This  drawing  was  transferred  to  a 
sheet  of  clear  plastic  and  inked  black  or  white  depending  on  the  predominant 
polarity  of  the  target.  The  transparency  was  photographed  against  a  neutral 
gray  background  at  a  distance  calculated  to  make  the  projected  target  the 
same  size  as  in  the  projection  of  the  complete  picture,  and  prepared  on  a 
70-mm  slide. 

The  targets  were  selected  at  random  in  respect  to  their  positions  within  the 
pictures  to  prevent  any  display -position  pattern  that  would  have  invalidated 
the  results  of  the  psychophysical  experiments.. 

SLIDE  PROJECTOR  AND  VIEWER 

As  shown  in  figure  1,  the  projector -viewer  had  three  12-  by  12-in.  screens 
arranged  side  by  side.  It  contained  two  70-mm  slide  projectors.  One  pro¬ 
jected  the  target  image  on  the  left-hand  screen  at  a  magnification  of  5.  5. 

The  other  projected  the  slide  of  the  entire  radar  scene  on  the  center  screen 
at  the  same  5.  5  magnification,  and  a  portion  of  the  same  slide  on  the  right- 
hand  screen  at  a  magnification  of  25;  this  was  done  through  two  channels  as 
shown  .n  figure  2. 

The  center  image  was  formed  at  a  set  ol  cross  hairs  on  a  translating  plate 
through  a  127 -mm  f/5,6  lens,  and  projected  onto  the  screen  through  a  180- 
m.m  f/4.5  lens.  The  right-hand  image  was  formed  at  another  set  of  cross 
hairs  on  the  same  translating  plate,  and  projected  onto  the  screen  through  a 
50 -mm  f/2.3  lens.  The  plate  was  driven  in  the  X  and  Y  directions  by  servo¬ 
mechanisms  controlled  through  an  ASQ-28  tracking  handle  operated  by  the 
subject.  When  the  plate  was  moved,  the  center  -  screen  radar  image  remained 
stationary  and  the  cross-hair  image  moved  over  it;  on  the  right-hand  screen, 
the  cross-hair  image  remained  stationary  and  the  radar  image  moved.  The 
center-screen  cross  hairs  were  used  by  the  subject  to  indicate  the  target;  the 
mag.Tified  image  on  the  right-hand  screen  simplified  the  checking  of  the  ac¬ 
curacy  of  the  target  identification.  The  cross-hair  images  were  inserted  in 
this  fashion  to  ensure  sharp  focus  and  reduce  parallax. 

The  projector -viewe  r  was  equipped  with  three  electrical  clocks  to  record  the 
times  consumed  in  studying  the  target,  locating  the  target  in  the  complete 
scene,  and  positioning  the  cross  hairs  over  the  target.  A  fourth  clock  measured 
the  total  time. 

TEST  SEQUENCE 

The  exj^erirnenter  seated  the  subject  at  the  projector -viewer  in  front  of  the 
three  screens.  The  experimenter  read  a  standard  set  of  instructions  to  the 
subject  (see  appendix  V).  Seven  practice  trials  were  run.  Then  the  experi¬ 
menter  projected  the  first  test  target  on  the  le'ft-hand  screen,  automatically 
starting  the  first  clock. 
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Figure  1  -  Slide  Projector- Viewer 


Figure  2  -  Schematic  of  Projector -Viewer  Center  and  Right-Hand  Channels 
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The  subject  studied  the  target  for  as  long  as  he  desired.  When  he  decided 
t.iat  he  had  memorized  it,  he  said  "Ready."  The  exper im.enter  then  projected 
che  complete  picture  on  the  center  screen,  an  action  that  simultaneously 
stopped  the  first  clock  and  started  the  second.  The  target  remained  on  the 

hand  screen.  The  subject  scanned  the  picture,  and  when  he  decided  that 
i  found  the  target,  he  pressed  a  button.  This  stopped  the  second  clock 
;.luminated  the  right-hand  screen.  The  original  target  and  complete  pic¬ 
ture  remained  on  their  respective  screens.  The  subject  positioned  the  center- 
screen  cross  hairs  over  the  target  he  had  identified  and  pressed  a  second 
button,  stopping  the  tliird  clock  and  completing  the  sequence.  The  accuracy 
c‘‘  the  identification  was  checked  by  means  of  the  magnification  on  the  right- 
hand  screen. 

The  experimenter  recorded  the  times  and  whether  the  target  had  been  identified 
correctly,  and  reset  the  clocks.  A  maximum  of  300  sec  was  allowed  for  locat¬ 
ing  the  target  (second  clock),  after  which  the  response  was  scored  as  incorrect. 
The  sequence  was  repeated  for  the  remaining  29  targets  and  scenes.  The  slides 
were  rearranged  in  a  random  sequence  between  each  test  series. 

RESPONSE  MEASUREMENTS 

At  the  time  of  these  human  response  experiments,  there  had  been  no  decision 
as  to  v/hat  measurements  would  be  used  in  formulating  the  prediction  metric 
(section  IV).  To  give  a  wide  choice,  the  complete  data  were  recorded  as  shown 
in  table  I.  As  it  turned  out,  only  the  average  recognition  time  for  each  slide 
was  used.  These  medians  are  tabulated  in  the  third  column  of  table  I,  and 
plotted  graphically  in  figure  3.  Plate  21  is  omitted  from  figure  3  because  none 
of  the  subjects  identified  this  target  within  the  300-sec  time  limit. 

Because  30  plates  and  30  subjects  were  used,  the  total  correct  and  incorrect 
identifications  would  be  expected  to  total  900.  They  total  only  867  because 
plate  8  was  found  to  be  unsatisfactory  after  seven  tests,  and  there  were  mal¬ 
functions  in  26  other  instances.  One  of  the  practice  plates  was  substituted  for 
the  defective  slide  and  designated  plate  8b. 

In  figure  4,  the  plates  are  arranged  in  groups  of  five  in  the  order  in  which  they 
were  viewed.  The  time  required  for  recognit  on  was  about  the  same  for  plates 
1  through  5,  11  through  15,  and  21  through  25,  and  about  the  same  for  plates  6 
through  10,  16  through  20,  and  26  through  30.  Thus,  there  was  no  indication 
that  learning  acquired  during  the  tests  affected  the  results. 
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TABLE  I 


ACTUAL  TARGET -IDENTIFICATION 
TIMES.  ORIGINAL  RADAR  PLATES 


Plate 

Identification  time  (sec) 

Number 

correct 

I 

i 

Percent 

incorrect 

T ime  to 
study  target 
(sec) 

Mean 

Median 

Standard 

deviation 

Percent 

correct 

Number 

incorrect 

19 

4  8 

5.  0 

1  4 

28 

93  3 

I 

2  1 

6  7 

6  9 

3 

5,  1 

9.  £  ; 

1  8 

30 

100.  0 

0 

0  0 

8.  6 

16 

6  6 

•:  0 

3.6 

28 

93.  3 

2  1 

6  7 

9  7 

7.2 

6.0 

3  9 

30 

100  0 

0  ' 

0  0 

9  3 

18 

7.7 

6.0 

6.0 

29 

96.  7 

1 

3  3 

8.  7 

4 

8  1 

6  0 

5.  7 

28 

96.  6 

1  ! 

3  4 

8.  7 

23 

8.  8 

6.0 

6.  7 

30 

100.  0 

0  : 

0  0 

8.  9 

14 

9.  8 

5.  0 

h  5 

30 

100  0 

0  ! 

0.  0 

7.  3 

7 

10.  1 

6.  0 

9.0 

29 

96.  7 

1 

3.  3 

10.  7 

17 

10.  8 

6  5 

9.2 

30 

100  0 

0 

0  0 

7.  3 

25 

12.  3 

12.  0 

7.7 

21 

84.  0 

4 

16  0 

9.9 

30 

14.7 

10.  0 

13.4 

27 

90.  0 

3 

10.  0 

12.  2 

15 

16.  8 

110 

14.2 

23 

92  0 

2 

8.  0 

9.  2 

20 

18.4 

11.5 

17  3 

20 

66.  7 

10 

33.  3 

10.  7 

24 

19.3 

7.5 

20.  4 

22 

75.9 

7 

24  1 

11  5 

27 

19  ? 

17.0 

15.  3 

18 

60.  0 

12 

40.  0 

6.  2 

5 

20.  7 

9.0 

29.  5 

28 

93.  3 

2 

6.  7 

7  2 

20.  8 

15  5 

19  7 

26 

86.7 

4 

13.  3 

9.  2 

28 

24.4 

14.  0 

33.  8 

27 

90.  0 

3 

10.  0 

9.  5 

26 

26  2 

21.0 

26  7 

25 

83.  3 

5 

16.  7 

8.  9 

13 

23  0 

30.  7 

27 

96.  4 

1 

3.6 

10.  5 

1 1 

34.  3 

26.0 

31. 8 

21 

72.  4 

8 

27.6 

1  ^ 

12 

37  9 

36  0 

20  9 

14 

60  9 

9 

39  1 

!  10. 6 

8b 

45  6 

34  5 

25  1 

16 

,  69  6 

7 

30.  4 

I  11.5 

6 

62  6 

52  5 

40  1 

24 

;  80.  0 

6 

20.  0 

7  3 

29 

64  1 

25.  0 

74  4 

12 

40  0 

18 

60  0 

7.  5 

1 

65  8 

49.  0 

55.6 

13 

48.  1 

14 

51  9 

9.  1 

2 

72.3 

47  0 

!  ■^7  3 

23 

1  76  7 

7 

23.  3 

8.  1 

22 

1  14 

145  0 

'  73  5 

1  1 

37. 9 

18 

62  1 

10.  0 

!  .  .  . 

0 

0. 

30 

100.  0 

11.3 

All 

1 

slides 

36.64 

21  46 

23.55 

i_ 

690 

78.  35 

J _ 

117 

21.65 

9  2 

i _ 
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MEAN  IDENTIFICa.TION  TIME  (SECONDS) 


Figure  3  -  Actual  Target-Identification  Times,  Original  Radar  Plates 


Figure  4  -  Order  of  Presentation  of  Original  Radar  Plates  versus 
Target -Identification  Times 
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SECTION  m 


SELECTION  AND  MEASUREMENT  OF  PICTURE  CHARACTERISTICS 


INTRODUCTION 

Duplicates  of  the  positive  plates  of  the  high  -  re  solution  radar  pictures  used  in 
the  human  response  experiments  were  used  to  establish  the  other  data  required 
for  the  prediction  metric  -  the  measurements  of  the  physical  characteristics 
that  were  believed  to  influence  the  target-recognition  times.  The  selection  of 
these  characteristics  was  based  on  many  considerations.  The  number  of  pos¬ 
sible  indices  of  the  characteristics  was  large,  and  it  was  reasonable  to  suppose 
that  many  would  have  low  correlation  with  the  human  response  data.  The  prob¬ 
lem  was  to  select  characteristics  that  were  believed  to  aifect  the  target-finding 
performance  and  that  could  be  measured  automatically  with  instruments  avail¬ 
able  under  the  contract.  A  further  limitation  required  that  such  qualities  as 
resolution,  emulsion  graininess,  and  acutance  be  ignored. 

The  detection  of  targets  in  high  -  resolution  radar  displays  presents  peculiar 
problems.  Ground-range  displays  were  used;  because  of  the  type  of  projection, 
these  have  more  resemblance  to  top>ographic  maps  than  to  optical  photographs. 
For  example,  objects  the  same  distance  apart  on  the  ground  appear  at  the 
same  relevant  distance  apart  on  such  radar  displays,  independent  of  how  far 
away  they  are  from  the  radar  sensor:  in  optical  photographs,  of  course,  the 
farther  the  objects  are  from  the  camera,  the  closer  together  they  appear  on  a 
picture.  Also,  the  brightness  of  an  object  illuminated  by  radar  is  considerably 
different  from  that  in  ordinary  photography.  Furthermore,  the  appearance  of 
an  object  changes,  often  radically,  with  even  small  changes  in  the  direction 
from  which  it  is  viewed  by  the  radar  because  the  reflection  angles  of  the  sur¬ 
faces  jf  the  objects  differ  with  the  direction  of  the  illuminating  beam. 

However,  certain  geometric  features  of  ground  objects  remain  the  same  in 
radar  pictures  regardless  of  changes  in  viewing  angle,  weather,  time  of  day, 
and,  to  a  great  extent,  the  season.  Such  features  are  valuable  guideposts  to 
which  an  observer  may  mentally  attach  local  geographical  coordinates.  In 
general,  they  have  strong  contour  sufficiently  distinguishable  by  the  sensor 
against  a  variety  of  backgrounds.  They  also  possess  connectivity  in  the 
psychological  sense,  which  preserves  their  identity  in  spite  of  noise  corruption 
and  gives  them  spatial  extent.  Runways,  land-w'ater  boundaries  such  as  shore 
lines,  rivers,  and  other  drainage  nets,  and  communication  lines  such  as  roads, 
railroads,  and  bridges  belong  in  this  category  along  with  metal  towers  of 
power  lines,  whose  truss  work  produces  nonchanging  radar  reflections,  towns 
and  cities  with  their  prominent  street  patterns,  and  industrial  complexes  and 
depots,  military  or  otherwise. 

In  spite  of  their  individuality,  even  these  targets  often  are  difficult  to  locate 
because  of  the  complexity  of  the  background  in  which  they  are  embedded. 

Their  visual  detection  and  identification  depends  eith'  r  on  characteristic  and 
outstanding  individual  invariants,  or  features  such  as  texture  that  may  be  ex¬ 
pressed  in  statistical  terms.  Even  invariant  features  may  not  make  a  target 
immediately  obvious  when  the  background  surrounding  it  contains  similar 
features  or  when  the  picture  scale  is  inappropriate.  Picture  interpretability 
under  these  circumstances  may  depend  upon  differences  betv/een  pertinent 
statistical  parameters  of  the  target,  on  the  one  hand,  and  of  the  background 
on  the  other. 
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SELECTION  OF  CHARACTERISTICS 


On  basis  of  the  above  considerations,  five  physical  characteristics  of  the  pic¬ 
tures  were  selected  for  measurement.  One  was  the  size  of  each  target.  The 
other  four  were  based  on  transmissivity;  that  is,  the  clearness  or  opaqueness 
of  the  positive  glass  plates  that  governed  the  passage  of  light  through  them. 
These  characteristics  were; 

The  average  transmissivity  of  each  target,  /i^,  in  relation 
to  the  average  transmissivity  of  the  area  immediately  sur¬ 
rounding  each  target,  //-p  (see  Measurement  of  Character¬ 
istics,  below) 

The  variance  (spread  about  the  mean)  of  the  transmissivity 
of  each  target,  s^,  in  relation  to  the  variance  (spread  about 
the  mean)  of  the  transmissivity  of  the  area  immediately 
surrounding  each  ta^'get,  s^ 

The  average  rate  of  change  (detail)  in  the  transmissivity  of 
each  target,  in  relation  to  the  average  rate  of  change 

(detail)  in  the  transmissivity  of  the  area  immediately  sur¬ 
rounding  each  target,  //q 

The  variance  (spread  about  the  mean)  of  the  rate  of  change 
(detail)  in  the  transmissivity  of  each  target,  s^j,  in  relation 
to  the  variance  (spread  about  the  mean)  of  the  rate  of  change 
in  the  transmissivity  of  the  area  immediately  surrounding 
each  target,  s^ 

These  were  dimensionless  quantities  expressed  in  ratios, 

MEASUREMENT  OF  CHARACTERISTICS 

The  transmissivity  of  each  radar  picture  was  measured  by  scanning  it  with  a 
flying  spot  scanner.  The  light  that  passed  through  the  positive  transparency 
was  focused  on  a  photomultiplier  tube  whose  output  was  a  voltage  that  varied 
in  direct  proportion  to  the  intensity  of  the  light.  This  voltage  was  fed  to  an 
anaiog  computer  programmed  to  process  it  mathematically  into  four  respec¬ 
tive  signals  representing  measurements  of  the  four  selected  transmissivity 
characteristics. 

A  flying  spot  scanner  was  built  for  the  measurements.  It  is  shown  in  figure  5. 
It  had  two  identical  channels,  one  to  scan  the  picture  and  generate  signals  to 
the  computer,  and  the  other  to  display  the  scan  to  the  scanner  operator  for 
monitoring  and  location  of  the  targets.  Generally,  in  these  devices  the  point 
of  light  (flying  spot)  is  generated  in  a  cathode-ray  tube  and  moved  across  the 
face  of  the  tube  in  paths  that  give  the  desired  scan  pattern,  or  raster.  Usually, 
this  scan  is  in  continuous  single  lines  traversing  horizontally,  as  in  TV  sets. 

The  high -resolution  radar  pictures  required  a  sequential  two -direction  ortho¬ 
gonal  scan  because  many  of  the  gradients  were  oriented  in  a  direction  other 
than  liorizontal.  This  sweep  voltage  was  programmed  by  the  computer  to  scan 
each  of  256  square  cells  of  the  picture  individually  10  times  and  10  times  ver¬ 
tically  in  the  sequence  shown  in  figure  6.  Figure  7  is  an  enlargement  of  the 
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Figure  6  -  Scan  Geometry 
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Figure  7  -  Scan  Pattern 
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scan  pattern  for  an  entire  picture.  The  single  dots  were  caused  by  crossings 
of  the  vertical  and  horizontal  scan  lines;  the  boundaries  of  the  256  cells  are 
marked  by  dots  caused  by  overlapping  sweeps.  The  distortion  was  due  to  the 
large  f-stop  of  the  lens  in  the  optical  path.  Figure  8  shows  how  the  flying  spot 
scanner  translated  an  optical  photograph  with  this  scan.  In  this  case,  the  light 
transmitted  through  the  positive  by  the  flying  spot  was  recorded  on  an  unexposed 
photographic  film. 

The  computer  generated  four  simultaneous  outputs  that  were  recorded  as  volt¬ 
age  variations  as  shown  in  figure  9.  Each  of  the  channels  represented  the 
average  values  of  one  of  the  transmissivity  characteristics  measured  cell  by 
cell.  A  fifth  channel  related  these  average  values  to  the  individvial  cells.  The 
measurements  were  dimensionless,  since  only  their  ratios  were  of  concern. 

The  scanner-computer  system  was  connected  and  programmed  to  perform  the 
following  four  operations  for  each  of  the  256  picture  cells: 


m 


j  =  1 

which  was  a  measure  of  the  average  (mean)  transmissivity; 


which  was  a  measure  of  the  variance  (spread  about  the  mean)  of  the  transmis¬ 
sivity; 


m 


which  was  a  measure  of  ihe  detail  (the  average  absolute  value  of  the  rate  of 
change  of  the  transmissivity);  and 
m 
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which  was  a  measure  of  the  variance  of  detail.  The  area  (a)  of  each  cell  oc¬ 
cupied  by  a  target  was  measured  by 


a.  = 
1 


x^)(y 


i  +  1 


(5) 


for  use  as  the  unit  of  measurement  for  target  size. 
In  the  above  equations, 

i  is  the  cell  number  in  the  cell  matrix; 


m  is  the  total  number  of  equidistant  scans  in  the  x  or  y  di¬ 
rection  within  each  cell  (in  this  case,  m  =  10); 


j  is  the  order  number  of  each  individual  scan; 


X.  ,  ,  -  X.  is  the  size  of  each  cell  in  the  x-axis  direction; 
1+1  1 


y.  ,  -  V  is  the  size  of  each  cell  in  the  direction  of  the  y 
axis  (the  two  sides  were  chosen  to  be  equal,  hence  x  = 


X  . 


i  +  1 


-  X.  =  y 


i  H-  1 


-  Yj  =  y); 


t(x,  y)  is  the  transmissivity  of  the  picture  depending  on  the 
position  X,  y  within  the  cell;  and 


3 1  3t 

- —  and  —  are  the  partial  derivatives  of  t(x,  y)  with  respect 
ax  dy  \  >  7  I 

to  X  and  y. 


On  the  basis  of  these  cell  values,  the  five  parameters  were  computed  sepa¬ 
rately  for  the  background  and  the  target.  For  example,  in  the  case  of  the 
average  transmissivity,  the  value  for  the  whole  picture  (the  sum  total  of  the 
N  =  256  cells)  was 


H  =  — 

T  N 

i  1 

and  for  the  target  alone  the  average  transmissivity  was 


n 


i  =  1 


(6) 


(7) 


In  the  expression  immediately  above,  only  the  n  cells  occupied  by  the  target 
were  considered.  The  other  transmissivity  parameters  were  determined  in 
similar  fashion: 


s 


ti  ’ 


(8) 
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1 


(9) 


and 


s 

t 


n 


n 


ti 


(10) 


(11) 


(12) 


(13) 


The  relative  size  of  each  target  to  the  entire  picture  area  (ail  256  cells)  was 
computed  on  the  basis  of 

=  N  icy 


and 


(14) 


=  n  xy  . 


Where  a  cell  was  filled  completely  by  a  portion  of  a  target,  the  computer  gave 
the  actual  average  of  each  of  the  four  characteristics  for  that  particular  cell. 
However,  the  periphery  of  each  target  cut  tnrough  many  cells  and  in  these 
cases  the  computer  did  not  distinguish  between  the  cell  target  area  and  the 
cell  background  area.  This  required  manual  analysis  and  computation  for 
each  such  cell  for  each  of  the  four  ci;aracteristics. 


Figure  10  is  an  illustrative  case  of  the  above.  Here,  is  the  portion  of  a 
cell  belonging  to  a  target,  and  1  -  X^  is  the  portion  of  the  same  cell  belonging 
to  the  background.  The  computer  found  the  average  transmissivity,  of  the 
entire  cell  by  means  of 


X  u  (1  -  X,)P  =  r, 

1  ti'  1  r  1 


with 
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Figure  10  -  Geometry  of  Target  Cell-Size  Measurement 

-  transmissivity  of  the  cell  target  area, 

-  transmissivity  of  the  cell  background  area,  and 

r.  =  the  computer  readout  of  the  average  for  the  entire 
^  cell. 

To  find  the  average  transmissivity  of  a  cell  target  area,  A  was  estimated  by 
visual  examination  and  th<i  following  was  solved  manually; 

r;  -  (1  -  A.)Mt 

= - V - • 

1 

Then  for  all  such  cells  (n). 


^  'p 
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Similar  manual  computations  were  required  for  all  such  cells  to  determine 
the  other  desired  ratios.  In  the  first  measurement  attempts,  positive  plates 
representing  18-  by  1 8 -nautical -mile  ground  areas  were  used.  Because  of  the 
distortion  limitations  of  the  flying  spot  scanner  optical  path,  these  pictures 
were  limited  to  2.25  by  2.25  in.  Thus  each  of  the  256  cells  into  which  each 
picture  was  divided  represented  a  ground  area  of  1.265  square  nautical  miles. 
At  this  scale,  the  cell  target  areas  (X  )  could  not  be  estimated  accurately. 

Also,  targets  of  the  airfield  category  Viad  to  be  excluded  because  it  \/as  all  but 
impossible  to  estimate  their  contributions  (for  this  and  other  reasons,  the 
original  30  plates  were  reduced  to  20).  Many  of  the  other  ratios  had  negative 
values,  as  shown  in  table  II,  that  indicated  unacceptable  errors  of  measure¬ 
ment  due  to  the  small  scale.  In  these  measurements,  the  background  trans¬ 
missivity  war  taken  as  the  average  for  all  256  cells.  This  average  was  com¬ 
puted  manually  from  the  recordings. 

A  new  scale  factor  was  adopted.  This  was  done  by  enlarging  each  target  and 
its  immediate  surrounding  area  six  times,  which  incidentally  was  about  the 
same  scale  as  the  projections  in  the  human  response  experiments.  The  20 
new  positive  plates  were  scanned  and  measured  as  before  by  the  flying  spot 
scanner  and  the  analog  computer,  and  the  required  manual  computations  were 
repeated.  At  this  new  scale,  each  of  the  2  56  cells  covered  a  ground  area  of 
only  0.  035  square  nautical  mile.  This  reduced  the  effect  of  target  and  back¬ 
ground  mixing,  because  now  the  targets  were  sufficiently  large  to  allow  better 
estimation  of  the  contribution  of  the  target  area  in  each  partially  occupied  cell, 
and  the  matrix  was  so  placed  that  there  were  fewer  cells  so  occupied.  How¬ 
ever,  even  at  this  scale,  airfields  had  to  bj  excluded  because  they  could  not 
be  measured. 

Two  sets  of  data  were  prepared  from  the  large-scale  pictures.  One  expressed 
the  ratios  of  the  four  transmissivity  characteristics  of  the  targets  to  the  back¬ 
ground  transmissivity  characteristics  averaged  over  all  256  cells,  as  shown 
in  table  III.  The  other  expressed  these  ratios  to  the  background  transmissivity 
averaged  over  the  immediate  background  of  each  target;  that  is,  those  cells 
partially  occupied  by  a  portion  of  a  target,  and  those  touching  the  edge  of  a 
target.  The  second  set  of  data  are  shown  in  table  IV.  Whether  a  limited  or 
overall  background  was  used  had  little  effect  on  the  transmissivity  ratios,  as 
shown  by  comparison  of  the  first  columns  in  the  tables. 

Unexpectedly,  it  was  found  that  the  four  transmissivity  characteristics  varied 
in  unison  in  almost  direct  proportion  from  cell  to  cell  in  many  of  the  pictures. 
This  is  shown  in  figure  11,  which  is  part  of  the  computer  recording  for  one  of 
the  large-scale  pictures.  The  scatter  diagram  for  versus  for  all  256 
cells  of  this  picture  (plate  17)  is  plotted  in  figure  12.  It  shows  a  high  degree 
of  curvilinear  correlation  of  the  tvpe 

''ti  =  ‘'"di  "  • 

where  k,  m,  and  c  are  constants.  Similar  correlations  are  indicated  by  fig¬ 
ures  13  and  14,  which  are  scatter  diagrams  for  versus  for  limited 

and  overall  target  backgrounds,  respectively. 

The  unexpected  correlation  between  these  characteristics  was  undesirable  be¬ 
cause  it  weakened  the  descriptions  of  the  pictures.  When  these  characteristics 
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were  selected,  they  had  been  expected  to  be  independent  of  each  other.  Time 
prevented  selection  and  measurement  of  new  characteristics.  Thus,  the 
measurements  discussed  above  were  used  in  deriving  the  human-response  pre¬ 
diction  equation  as  discussed  in  section  IV. 

TABLE  U 

CHARACTERISTICS  OF  TARGET  VERSUS  OVERALL 


BACKGROUND,  ORIGINAL  RADAR  PLATES  AT  SMALL  SCALE 


Cell  location 

Ratios  of  transmissivity 
characteristics^ 

Plate 

Column 

Row 

•t  't 

*d  *D 

1 

6 

1 1  (0.  35),  12  (0.  10) 

-0.  5 

0.6 

3.3 

2 

\ 

13  (0.  45) 

.  .  . 

a  a  a 

a  .  a 

3 

1  1 

8  (0.  80) 

1.6 

1.5 

1.9 

4 

9 

3  (0.  10) 

2.3 

31.9 

3.9 

18.7 

5 

6 

9  (0.  65),  10  (0.  05) 

.  .  . 

a  a  a 

a  a  a 

a  a  a 

6 

8 

5  (0.25) 

.  .  . 

a  a  a 

a  a  a 

a  a  a 

7 

7 

8  (1. 00) 

1.4 

1.  1 

1.4 

1.6 

8 

7 

6  (0.25).  9  (0.  35) 

1.5 

0,2 

1.9 

2.  0 

9 

5 

4  (0.  50) 

.  .  . 

a  a  a 

a  a  a 

a  a  a 

10 

7 

5  (0.  40) 

1.7 

1.3 

1.7 

1  1 

9 

2  (0.  10) 

0.  8 

-0.4 

-1.0 

‘2 

7 

9  (0.  90) 

1.  1 

1.2 

1.2 

13 

6 

12  (0.  03) 

8.  8 

12.  3 

13.9 

14 

7 

10  (0.  15) 

.  .  . 

a  a  a 

a  a  a 

15 

12 

2  (0.  50) 

1.2 

1.3 

2.  1 

16 

6 

9  (0.  50) 

1.3 

1.9 

2.  1 

17 

8  (0.  05)* 

9 

a  •  a 

a  a  a 

a  a  a 

a  a  a 

9  (0.70) 

.  .  . 

.  .  . 

.  .  . 

.  .  . 

.  .  . 

18 

2 

4  (0.  15) 

-0.  1 

4.  1 

-0.  0 

-0.2 

19 

4 

4  (0.  20).  5  (0.  10) 

a  .  a 

a  a  a 

a  a  a 

a  a  a 

20 

5 

4  (0.  10) 

a  .  a 

a  a  a 

a  a  a 

a  a  a 

8 

7  (0.  10).  8  (0.  50) 

.  a  a 

a  .  a 

a  a  a 

a  a  a 

21 

9 

7  (0.  10) 

1.7 

2.8 

2.7 

3.2 

12 

10  (0.  005) 

a  a  . 

a  a  a 

a  a  a 

a  a  a 

22 

1 1 

5  (0.  10) 

1.8 

18.9 

2.  3 

3.  3 

23 

1 

15  (0.  10) 

-1.7 

-7,  7 

-2.2 

-2.  0 

24 

6 

8  (0.  07) 

a  a  a 

a  a  a 

a  a  a 

a  a  . 

7 

9  (0.  10) 

a  .  . 

a  a  a 

a  a  a 

a  a  a 

25 

3 

13  (0.  20) 

0.  0 

2.9 

-0.  5 

-0.  7 

26 

9 

7  (0.  08) 

-3.  1 

8.  5 

-3.  3 

-2.  7 

27 

5 

5  (7.  18) 

1.5 

13.  3 

2.  2 

2.4 

28 

7  (1.0) 

7 

1.6 

1.4 

1.6 

1.7 

8 

.  .  . 

a  a  a 

a  a  a 

a  a  a 

a  a  a 

29 

8 

8  (0.  05) 

a  a  a 

a  a  a 

a  a  a 

a  a  a 

30 

8 

6  (0.  50).  7  (0.  10) 

1.2 

-3.  2 

-3.  1 

_ 

1.2 

Numbers  in  parentheses  show  percentage  o{  cell  occupied  by  target. 


*  target  to  background;  s^  s.^  *  target  variance  about  mean  to  background  variance 
about  mean;  ~  target  rate  ot  change  (detail)  to  background  rate  of  change,  * 

target  detail  variance  to  background  detail  variance. 


19 


TABLE  m 


CHARACTERISTICS  OF  TARGET  VERSUS  OVERALL 
BACKGROUND,  ORIGINAL  RADAR  PLATES  AT  LARGE  SCALE 


Hlale 

Ratios  of  t  ranBiiii  s  s  1 V  it  y 
c  li.i  rai  t «■  r  1  St  1 1  s 
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ba i  k r o ii  1 1  ci  i  h a  r a i  1 1-  r  i  .s t if 

rgrl  a 
lo 

H 

i)(i 

Target  si/.f 

(a  X  256) 

Mt 

a 

1 

n 

“,1 

s^  -  s  1, 

DBI 

Bl 

-  "o 

M,. 

s,. 

B 

"l> 

UH 

s  ,, 

s 

n 

1 

1  S3 

1  42 

1  t>5 

1  .  30 

0  53 

0  4  2 

0  65 

0. 

3  7 

a 

.  .  . 

.  . 

.  .  . 

.  .  . 

.  .  . 

.  .  . 

.  .  . 

• 

.  .  . 

3 

1  99 

2  30 

1  4  1 

1  44 

0  00 

1  30 

0.4  1 

0 

44 

5  1 

4 

3  00 

3  02 

1  7  1 

1  04 

2  00 

2  02 

0.7  1 

0 

04 

1  0 

') 

6 

: : 

: : : 

:  : 

:  :  : 

'  :  : 

:  :  ‘ 

: 

; ; 

7 

0  68 

4  63 

0  06 

1  22 

0  32 

3.63 

0  04 

0 

22 

7  .  7 

8 

1.01 

1.10 

1  47 

1  45 

0  01 

0  1 0 

0.  47 

0. 

45 

8  0 

9 

•  .  . 

«  •  . 

.  .  « 

•  .  • 

•  •  • 

•  •  » 

.  .  • 

* 

.  . 

.  .  . 

1  0 

1  1  5 

0.  80 

0.  94 

1  .  C  1 

0.  15 

0.  20 

0  06 

0. 

01 

3.  3 

1  1 

1  41 

I  .  39 

1.45 

2  05 

0  41 

0  30 

0  4  5 

1 

05 

2  0 

U 

1  35 

1  .  52 

1  3  1 

1  33 

0  35 

0.  52 

0  31 

0 

33 

1  3  5 

1  3 

1  77 

3  00 

1 . 07 

1  29 

0  77 

2  00 

0.  07 

0 

29 

0  8 

14 

.  .  . 

.  .  . 

.  .  . 

.  .  . 

.  .  * 

*  .  . 

.  .  . 

. 

»  . 

•  .  . 

1  5 

2  22 

2  50 

1  54 

1  63 

1  22 

1  50 

0.  54 

0 

6  3 

3.  6 

16 

2.  10 

9  6 

1  44 

2  51 

1  10 

8.  6 

0  44 

1 

51 

5  0 

17 

.  .  . 

.  .  . 

.  .  . 

.  . 

.  .  . 

.  .  . 

.  .  • 

. 

.  . 

.  .  . 

1  8 

141 

2  85 

0.  7  0 

1  22 

0  41 

'  .  85 

0  21 

0 

22 

1  5 

19 

.  .  . 

.  .  . 

.  .  . 

,  .  . 

.  .  . 

.  .  . 

.  ,  . 

20 

1 . 08 

1  99 

0  05 

1  14 

0  08 

0  00 

0  0  5 

0 

14 

1  5  7 

21 

•  •  ■ 

.  .  . 

.  .  . 

.  ,  . 

.  .  . 

.  .  . 

.  .  . 

. 

•  . 

22 

0  87 

0.  58 

0  73 

0  54 

0  13 

0  42 

0  27 

0 

46 

0.  6 

23 

3. 64 

7  72 

2  22 

2. 94 

2. 64 

6  72 

1  .  22. 

1 

04 

1  0 

24 

•  •  . 

.  .  • 

•  .  . 

*  •  • 

•  .  , 

.  .  • 

.  .  . 

. 

.  . 

25 

3  07 

2.  83 

5  54 

8.  22 

2  07 

1  83 

4  54 

i 

22 

1  ti 

26 

1  .  50 

2  82 

1 . 09 

1  60 

0  SO 

1  82 

0  00 

0. 

60 

0  0 

27 

1 . 43 

6  38 

1 . 86 

2  1  5 

0  43 

5  38 

0  86 

1 . 

1  5 

2  5 

28 

1 . 46 

1.6- 

1  22 

1  22 

0  46 

0.  68 

0  22 

0. 

22 

24  4 

29 

.  •  . 

.  .  . 

.  .  . 

•  .  . 

.  ♦  » 

.  . 

.  .  . 

. 

.  .  . 

30 

1  16 

4  00 

1  49 

1 . 64 

0  16 

3.  00 

0.  49 

0 

64 

12  8 
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TABLE  IV 


CHARACTERISTICS  OF  TARGET  VERSUS  IMMEDIATE 
BACKGROUND,  ORIGINAL  RADAR  PLATES  AT  LARGE  SCALE 


I 

Ratios  ol  transmissivity 
i'  ha  rat  t  e  r  1  st  u  s 

Ratios  of  diffor«‘iu  t-  between  tarj>et  and 
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“t 
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%1 

1 
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1 

1 
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1 

1 
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1 

II 

Bl 

■  "d 

1 
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(o  X  2  5tj) 

"l 

S.J, 

**13 

1 

r 

1 

1 

1 

1 

"d 

1 

1 

S 

D 

1 

1 

1  .  34 

1  .  1  S 

1  57 

1  1  1 

0  34 

0  1  5 

0  57 

0 

1  1 

3  7 

z 

.  .  . 

.  .  . 

a  a  a 

a  a  a 

a  a  a 

...  j 

a 

3 

2  04 

2  50 

1 . 44 

1  50 

1 . 04 

1  bO 

0  44 

0 

50 

5  1 

4 

5 

6 
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2  82 
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1  85 

1 . 98 

1  82 

0  66 

0 

85 

1  t 

•  •  a 

a  a  a 

: 

: 

: 

:  ’  : 

; 

7 

0.68 

4.  56 

0  98 

1  23 

0  32 

5  56 

0  02 

0 

23 

7  7 

H 

1.01 

1  1  1 

1  25 

1  19 

0  01 

0  1  1 

0  25 

0 

1  u 

8  0 

9 

.  .  . 

a  .  a 

a  a  a 
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a  a  . 

1  0 
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0  85 

0  91 
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0  16 

0  15 

0. 

09 

3  3 

1  1 

1  48 

1  48 

1  43 

2.01 

0  48 

0  48 

0.43 

1 

01 

i.  0 

12 

1  .  1  S 
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1  1  3 

0  89 

0.  IS 

0.  54 
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0 
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1  3 

1  .  S7 

2. 47 

0  89 
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0 

04 

0  8 

H 
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a  > 

.  a 
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0  57 

0. 
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16 
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1 

32 
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0 
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27 

1  .  24 
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0  55 

0 
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28 

1  .  38 
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F igure  1 1 


-  Analog  Computer  Recordings  Showing  Correlation  between 
Cell  Parameters  (Original  Radar  Plate  17) 
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(MILLIMETERS 


TARGET  DETAIL  BACKGROUND  DETAIL 


TARGET  TRANSMISSIVITY  BACKGROUND  TRANSMISSIVITY  i-j) 


Figure  14  -  Scatter  Diagram  of  Transmissivity  and  Detail  Ratios,  Overall 
Background,  Original  Radar  Plates  at  Barge  Scale 
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SECTION  IV 


SELECTION  OF  THE  HUMAN  PERFORMANCE  PREDICTION  EQUATION 


INTRODUCTION 

With  the  human  response  data  established  and  the  physical  characteristics  of 
the  pictures  measured,  the  problem  was  to  define  a  metric  that  would  relate 
these  two  sets  of  data  in  a  prediction  equation.  Three  metrics  were  defined 
by  three  teams  working  independently.  Metric  I  related  the  two  sets  of  data 
in  a  logarithmic  equation  as  discussed  below.  Metrics  II  and  III  were  predi¬ 
cated  on  the  same  general  type  of  picture  characteristics  as  used  in  Metric  I 
but  the  relationships  were  utilized  differently,  as  explained  in  appendixes  I 
and  II,  respectively.  Lack  of  time  prevented  attempts  to  validate  Metrics  II 
and  III. 


SYNTHESIS  OF  METRIC  I 

In  synthesizing  Metric  I,  the  intention  was  to  relate  the  four  transmissivity 
characteristics  of  target  and  background,  as  well  as  target  relative  size,  to 
the  psychophysical  scores  of  the  human  observers.  The  following  correspond¬ 
ing  characteristics  for  target  and  background  in  the  form  of  ratios  were  chosen: 

L  L  lA  and  . 

''T  ®T ’  ^d’  ®D  ’  “t 


The  aim  was  to  combine  these  ratios  in  a  manner  to  bring  out  the  differences 
between  target  and  background  in  both  brightness  and  detail.  Two  such  ways, 
both  attractive  for  their  computational  simplicity,  were  considered.  The 
first  one  arranged  each  parameter  pair  (such  as  and  in  the  form  of  the 

sum  of  their  ratio  and  the  inverse  of  the  ratio.  With  the  general  terms  b  and 
B  denoting  any  one  of  the  parameter  pairs,  then 


c  . 


This  arrangement  was  required  because  an  observer  attempting  to  distinguish 
the  target  from  its  background  would  benefit  as  much  from  b  being  greater 
than  B  as  from  the  reverse. 

For  the  cases  of  the  four  transmissivity  characteristics,  the  value  b  for  the 
target  and  the  value  B  for  the  background  could  be  such  that 

b  I  B  . 

This,  of  course,  would  not  be  true  about  the  target  size,  which  would  be  sub¬ 
ject  to  the  condition 


b  <  B  . 

The  purpose  of  the  constant  c  was  to  establish  an  appropriate  scale  origin  for 
measuring  /3.  On  this  basis,  the  following  five  dimensionless  parameters  were 
established: 


Z5 


T 


s 

T 


c  , 


Ilf.  _ 

'^5  —  - +  - +  c  , 


®d  ®D 

s  =  - +  - +  c,  and 

*  =d 


a  = 


ot 

^  4- 

— —  +  c 

o 


It  was  assumed  that  the  contribution  of  each  of  the  transmissivity  parameters 
would  become  zero  when  target  and  background  were  the  same  as  far  as  a 
particular  parameter  was  conf.erned;  thus 

/3  =  ^+  -^+c  =  0 


when 


_b 

B 


1  . 


From  these  two  equations, 


c  =  -2  . 

For  che  area  ratio,  however,  a  had  to  be  zero  when  the  s '  ze  of  the  target 
was  zero.  This  made 


c  =  0  . 
o 

When  these  values  of  c  and  c^  were  substituted  in  the  above  equations,  then 


„  (1  -  ^/V‘ 

r  V'-T 


(1  -  St/  S^)' 

t^®T 


s/s 


= 


(1  - 
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and 


together  with  the  condition  that  if 


then 


and 


because  none  of  the  dimensionless  parameters  would  have  any  meaning  when 
the  target  area  was  reduced  to  zero. 

The  other  way  of  arriving  at  similar  results  was  by  defining 


which  was  equivalent  to 


and,  therefore,  to  saying  that  was  the  ratio  of  the  absolute  value  of  the  dif¬ 
ference  in  transmissivity  characteristics  between  target  and  background  to 
the  value  of  the  same  characteristic  for  the  background.  Thus, 
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‘5  = 


1  - 


■ 


1  - 


D 


and 


Under  the  existing  state  of  knowledge  on  perception  and  recognition,  there  was 
no  psychophysical  law  that  could  dictate  how  the  physircal  parameters  discussed 
above  affect  human  operator  performance  in  target  detection  tasks.  However, 
the  values  of  the  purely  physical  characteristics  represented  by  the  parameters 
had  to  be  brought  into  some  quantitative  relationship  with  the  scores  obtained  in 
the  psychophysical  experiments. 

It  was  assumed  that  a  group  of  observers,  subjected  to  identical  briefing  and 
put  into  an  identical  experimental  environment,  when  presented  with  the  same 
picture  containing  a  particular  target  would  take  different  lengths  of  time  to 
identify  the  target  and  that  some  of  them  would  identify  the  wrong  target  some¬ 
times.  The  two  outcomes  of  primary  interest  in  such  experiments  are  the  per¬ 
cent  of  correct  decisions,  p  ,  and  the  time  to  reach  a  correct  decision,  s^. 
Only  the  latter  was  retained  for  correlation  with  the  physical  picture  charac¬ 
teristics  already  discussed.  It  was  assumed  that  this  correlation  could  be 
expressed  mathematically  in  the  form 

^  ,  sg) . 

c 


Because  available  knowledge  did  not  suggest  a  particular  form  of  the  function 
f(/x^,  s^,  ,  sg  ),  a  heuristic  approach  was  employed.  Forms  of  the  equation 

tried  wereTinear,  vectorial,  logarithmic,  and  exponential  in  the  four  trans¬ 
missivity  parameters,  while  the  relative  target  size  was  either  included  or 
excluded  in  each  case.  Specifically,  the  following  equations  were  tried: 


(A) 

(B) 


— ^  =  k,  +  k-,s  +  k-jfi  +  k.s  +  kc  , 

s^  1  r  z  T  3  5  4  ^  5 


=  a(k^M^  +  ^2®  7"  ^  ^3  ^5  ^  ^4  ^  ^  ^5  ’ 

c 


(15) 

(16) 


(C) 

(D) 


4®5  ^5  ’ 


r:  =  kj , 


(IV) 

(18) 
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=  k  j/n  k^ins^  +  ^  ^5  ' 

— ^  =  Q;(k,/n/^  +  k^fns  +  k^/n  +k.iVis-)  +  k-.  , 

s^  '  I  T  2.  r  354s'  5 

^  =  exp(k,^^+  k^s^  +  ‘'3''5+  *^4^5  ^  ‘'5)  • 

* 

—  =  exp  Cl(^y  ^  ^3  ^5^  ^4^5)  ^  ^3 


where  the  ratios  had  the  form  of 


(1  -  h/B)‘ 

b/B 


(19) 

(20) 

(21) 

(22) 

(24) 


or 

(b)  =  I*  'll'  >“> 

Hence,  the  total  of  alternative  expressions  was  16:  (A^),  (A^),  (B^),  (B|j^), 

(C^).  (C^),  (D^),  (D^),  (E^),  (E^).  (F^),  (F^),  G(a).  (G^,),  (H^),  and  (H^). 

Two  points  should  be  explained  at  this  juncture.  First,  ^/s^,  was  used  in  all 
these  expressions  and  not  s^  because  the  target  recognition  time,  s^,  was  ex¬ 
pected  to  increase  as  the  picture  parameters  /3  decreased.  This  suggested  a 
type  of  inverse  relationship  between  s^,  and  the  /3's, 

Second,  whenever  the  target  size  was  included  in  the  metric,  it  was  used  not 
as  an  additive  term  but  as  a  multiplier.  This  was  done  because  of  the  special 
role  target  size  played  here:  if  the  target  area  diminished,  recognition 
eventvially  v/ould  become  impossible  regardless  of  the  values  of  the  other 
parameters.  Also,  regardless  of  the  value  of  the  ratio  a,  no  recognition 
would  be  possible  if  all  the  other  parameters  dealing  with  transmissivity 
characteristics  were  zero,  because  in  this  case,  no  structural  feature  would 
be  present  to  distinguish  the  target  from  each  background. 

Consequently,  a  functioned  as  a  multiplicative  factor  differing  from  s^,  , 

and  Sg,  which  could  be  taken  additively.  There  was  evidence  (references  2 
through  5)  that  recognition  time,  s^,  increases  with  display  size,  Q-p,  and 
decreases  with  target  size,  a^.  The  reverse  behavior  had  been  noted  for  the 
probability,  p^,  of  correct  decisions  (reference  6). 

A  reasonable  approach  toward  ascertaining  the  validity  of  the  above  postulated 
functional  relationships  between  the  picture  physical  parameters  and  the  psy¬ 
chophysical  scores  was  the  use  of  multiple  regression,  a  purely  empirical 
procedure  for  deriving  a  mathematical  relationship  between  the  quantities  of 
interest  (/3's  and  s^).  It  is  often  useful  because  it  yields  estimates  and  predic¬ 
tions.  It  should  be  used,  however,  with  the  clear  understanding  that  the 
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relationship  is  a  purely  mathematical  devi-'e  devoid  of  any  specific  cause- 
and-effect  implications.  The  fact  that  a  number  of  variables  tend  to  vary 
together  according  to  some  mathematically  desirable  regularity  does  not  in 
itself  define  the  specifics  of  the  mechanism  involved.  Multiple  regression 
of  the  psychological  variable  s^.  on  the  picture  parameters  ,  s  ,  ,  s 

was  the  goal.  It  required  the  application  of  the  least -squares  method  to  the 
expression 


TVT  N 

=  1  7  _  , 


1=1*-  J 


^  'k,  +  k_s  +  k, /I  +  k.s 

s  V  1  r  3  4  5 

c.  \  Oj 


for  the  A,  case.  Similar  expressions  were  used  for  the  rest  of  the  cases, 
through 


Here,  ®7-j»  gj>  Sgj  were  the  values  of  the  transmissivity  characteristics 

for  the  j-th  plate  and  s^j  was  the  average  score  of  target  recognition  for  all 
observers  for  the  same  j-th  plate.  N  was  the  total  number  of  plates.  The 
regression  computations  were  performed  with  the  right-hand  member  of  the 
tried  metric  equations  A  through  H  always  in  a  linear  form  with  respect  to 
the  constants  kp  k,,  k^#  k^,  and  k^.  This  meant  that  in  the  C  case  the  met¬ 
ric  was  rearrangecf  as 


-tj  =  +  k_,Sg^  t  kg 

s 

c 

2 

before  being  introduced  in  the  expression  for  the  total  error  squared,  f  . 
Also,  in  the  G  case  the  metric  was  rearranged  as 


c 


Similar  rearrangements  were  applied  to  the  D  and  H  cases. 

The  least- squares  technique  yielded  for  a  system  of  five  equations  with  five 
unknowns,  kj  through  k^.  Solution  of  the  system  provided  the  numerical  values 
of  these  unknowns.  The  system  was  obtained  by  equating  to  zero  the  derivatives 
of  with  respect  to  each  of  the  unknowns,  kj  through  k^;  that  is. 


a 


air 

1 


0 


with 


i  =  1,  2,  3,  4,  5  . 

All  operations  involved  in  the  multiple  regression  were  performed  by  a  digital 
computer.  The  criterion  used  to  measure  the  goodness  of  fit  in  each  case  was 
the  correlation  coefficient,  r,  between  target  recognition  times  observed,  s^,q, 
and  the  corresponding  times  predicted,  by  the  equation.  By  definition, 

the  correlation  coefficient  was 


30 


(26) 


R  = 


N 


E 

i  =  1 


(s  -  s  ).{s  -  s  ). 

CO  CO  j  ce  ce  j 


no  a 

ej 


where  s  and  a  ,  and  s  and  a  were  the  average  and  standard  deviations  of 
CO  o  ce  e 

s  and  s  ,  respectively,  and  N  was  the  number  of  correlated  scores.  The 
CO  ce 

values  of  R  were  computed  as  shown  in  table  V. 


The  multiple  -  regres sion  technique  by  the  least-squares  method  showed  that 
Alternative  (equations  19  and  24)  most  nearly  fitted  the  target-finding  times 
recorded  in  the  original  human  response  experiments.  It  had  an  R  of  0.  689. 
Thus,  the  human  performance  prediction  equation,  with  the  five  k's  replaced 
by  their  numerical  values,  became 


■^  =  0.  00480  +  0*  0041  1  log^  s^  -  0.  00476  ' 

c 

0.00954  logg  Sg  +  0.05993  sec"^  .  (27) 


Equation  27  was  applied  individually  to  the  20  original  human -response  pic¬ 
tures;  the  results  are  plotted  in  figure  15  and  listed  in  table  VI. 

The  hypothesis  that  in  the  population  the  m.ultiple  correlation  was  zero  (that 
is,  that  the  obtained  R  was  not  significantly  different  from  0)  was  tested  by 
the  F  distribution  with 


_  R^  N-M-1 

^M,  N-M-1  "  J  _  M  ’ 

where  M  was  the  number  of  predictor  variables  (M  =  4).  The  value  of 

F.  .  X  .  1  computed  with  this  formula  was  F  .  ,  ^  =  3.  37.  On  the  other 

M,  N-M-1  ^  4,  15 

hand,  F^  at  the  0.05  level  of  significance  was  3.06;  hence,  the  hypothesis 

that  the  population  multiple  correlation  was  zero  was  rejected  at  the  0.  05  con¬ 
fidence  level. 
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TABLE  V 


COMPUTER  REGRESSION  ANALYSIS  OF 
ALTERNATIVES  AND  CONSTANTS  TO  ESTABLISH  METRIC  I 


Alt*-  rna  - 

k 

k  , 

k.. 

k  c 

R 

t  1  V  f 

1 

1 

3 

4 

8 
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ba  C  round 

A 
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t  0. 00b 3  3 
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-0.  0030 1 
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0.  380b 

B 
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C 
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i) 

t  0.  000^(y 
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X 

10' 

s 
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-0.  008  33 

f  0.  0083  3 

0.  7446 

t  0.  00480 

t  0.  004 1  1 

-0. 004 7b 

-0  00084 

♦  0.  0800  3 

0.  688*1 

a 

(1 

t  0.  0044Z 

*  0.  00  3  1  0 

-0  00874 

-0  00307 

♦  0.  074  3  8 

-0. 0478 

♦0  S4701 
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(1 

a 
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^b 
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TABLE  VI 


METRIC  I  APPLIED  TO  ORIGINAL 
RADAR  PLATES  TO  VERIFY 

regression  analysis 


Time  to  identify  target 
(sec ) 

Plate 

Actual 

Predicted 

1 

65.  8 

13.  3 

3 

5.  1 

11.8 

4 

8.  1 

13.  2 

7 

10.  1 

8.  6 

8 

45.6 

41.4 

10 

7.2 

7.8 

1 1 

34.  3 

17.9 

12 

37.9 

12.  9 

13 

30.4 

17. 0 

15 

16.  8 

12.  7 

16 

6.6 

13.  1 

18 

7.7 

9.7 

20 

18.  4 

12.  8 

22 

11.4 

42.  3 

23 

8.  8 

14.  3 

25 

12.  3 

26.  1 

26 

26.  2 

11.2 

27 

19.7 

15.  0 

28 

24.  4 

10.  9 

30 

14.7 

15.  0 
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SECTION  V 


GENERALIZATION  OF  METRIC  I 


RADAR 

Twenty  radar  photographs  were  selected  to  check  the  validity  of  equation  27. 
Thirteen  were  new  ones  selected  from  Quick  Check  flights;  the  other  seven 
were  carried  over  from  the  original  human  response  experiments  for  direct 
comparison  of  the  performance  of  two  groups  of  subjects. 

Positive  glass  plates  were  prepared  of  the  13  new  pictures.  The  same  charac¬ 
teristics,  except  for  target  relative  size  which  was  not  a  factor  in  the  human- 
response  prediction  equation,  were  measured  as  before  by  the  flying  spot 
scanner -analog  computer  system  and  manual  calculations. 

The  characteristics  of  three  of  the  new  radar  plates  (Nos.  3,  11,  and  13)  had 

negative  values  or  could  not  be  measured  at  all.  These  plates  were  discarded, 
reducina  tht;  number  to  10,  in  addition  to  the  seven  held  over  from  the  original 
human-response  experiments. 

Table  VTI  li  :is  '.he  numerical  values  Oj.  the  characteristics  and  the  predicted 
target- recognition  times  obtained  by  substituting  these  values  into  equation  27. 

At  the  sp'no  time  that  the  predicted  target-identification  times  were  being  es- 
tablisi  .^i  ^  ^  the  new  plates,  actual  recognition  times  were  being  recorded  for 

duplic;=' ■ .  plates  in  experiments  with  22  SAC  radar  navigators  at  Lockbourne 
Air  Force  Base,  Ohio,  The  equipment  and  procedures  were  the  same  as  in 
the  original  experiments  at  Griffiss  and  Pease  AFBs  except  that  only  five 
practice  trials  were  run  instead  of  seven.  These  five  practice  slides  were 
selected  from  those  used  at  Griffiss  and  Pease.  At  Lockbourne,  all  13  of 
the  new  high  -  resolution  pictures  were  used  because  it  was  not  known  at  that 
time  that  the  characteristics  of  three  of  them  could  not  be  measured.  These 
experiments  also  included,  of  course,  the  seven  plates  held  over  from  the 
original  tests. 

Table  VIII  shows  the  results  of  the  Lockbourne  experiments.  These  results 
are  plotted  in  rank  order  in  figure  16.  As  in  the  original  experiments  and  for 
the  same  reasons,  no  data  were  collected  in  some  tests. 

Figure  17  shows  the  correlation  of  the  identification  times  between  the  original 
and  validation  subjects  for  the  seven  slides  common  to  both  tests.  The  cor¬ 
relation  was  0.97;  it  was  significant  at  the  0.01  confidence  level. 

Table  IX  directly  compares  the  predicted  and  actual  target  identifications  for 
the  10  new  high-resolution  radar  plates.  The  correlation  coefficient,  R,  was 
0.  284,  which  was  not  significantly  different  from  zero  for  this  number  of 
plates.  When  plate  number  two  was  omitted,  the  correlation  betwetn  the  ob¬ 
tained  v'  ean  scores  and  the  predicted  scores  for  the  remaining  nine  plates 
was  only  -0.  034.  When  the  median  scores  were  used  instead  of  mean  scores, 
the  correlation  for  the  nine  plates  was  -0.048. 

The  difficulty  in  measuring  the  target  cells  with  adequate  precision  took  such 
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TABLE  VU 


METRIC  I  APPLIED  TO  NEW  RADAR  PLATES  TO 
PREDICT  TARGET -IDENTIFICATION  TIMES 


Ratlob  of  tranarniHBivity 
c  ha  ra  (.  t «'  ri  8t u  h 

Ratios  of  (J  if  f  e  ri'Mce  N 
h«*twf<*n  targ*‘l  and  b.uk- 
groiifid  1  h.irat  It'ristii  a  to 
bai  kg  round  i  lia  rac  t  o  r  i  t»t  ic  h 

Pr«‘du  ted 
ident  i  f  11  a 

■‘l 

“d 

Plate 

"t- 

8 

r 

tion  timr 
(sec  ) 

1 

1 . 44 

0  6S 

1  20 

0  70 

0  1399 

0  1886 

0.  0333 

0.  1286 

12.6 

1 

1  73 

2. 42 

1 . 49 

1  72 

0  3 1 09 

0.  8332 

0.  1611 

0  3014 

13  6 

4 

0.  67 

o  • 

. 

0  47 

0  63 

0  1626 

0.  3467 

0.  6977 

0  2173 

16  7 

5 

1 . 09 

1  .  68 

1.  03 

1  08 

0.  0074 

0  2129 

0  009 

0  0069 

8.  9 

6 

0.  28 

2.  98 

0  60 

0.  67 

1  8614 

1  3166 

0.  6000 

0.  3244 

12.8 

7 

1.31 

0.7  3 

0.76 

0  79 

0.  07  34 

0  0999 

0  0833 

0.  0668 

12  7 

8 

2  26 

1  .  49 

1  61 

1  46 

0  7026 

0  1611 

0.  1723 

0  1  397 

12.8 

2  71 

2.  30 

2.  06 

2  26 

1 . 0790 

0  7348 

0.  6378 

0.  6944 

16  4 

10 

123 

1  76 

1 . 01 

113 

0  04  30 

0  3228 

0  0001 

0.  0160 

9.  8 

1  1 

12 

1  3 

2.7  3 

8.  84 

[.  .  . 

1.  66 

2  14 

0  096  3 

6.96  3 

0  2010 

0.  6073 

12.4 

TABLE  Vm 

ACTUAL  TARGET -IDENTIFICATION  TIMES, 
NEW  SET  OF  RADAR  PLATES 


Identification  recognition 
I ime  ( »cc  ) 


Plate 

Mean 

1  Median 

[  Standard  j 
1  deviation 

Number 
cor  ret  1 

Percent 

correct 

I 

Number 

incorrect 

Here  ent 
incur  rec  t 

8 

7 

.  8 

6 

6 

6 

8 

22 

1  00 

0 

0 

0 

12 

8 

6 

1  4 

0 

1  1  . 

6 

22 

1  00 

0 

0 

0 

19 

9 

4 

1 

0 

4 

21 

1  00 

0 

0 

0 

14 

10 

0 

b 

0 

8. 

3 

22 

1  00. 

0 

0 

10 

10. 

2 

1  6 

•- 

9. 

7 

22 

100. 

0 

0 

6 

12 

1 

8. 

0 

10. 

0 

21 

96. 

6 

1 

4 

6 

17 

14 

9 

9 

0 

16 

6 

22 

1  00 

0 

0 

0. 

0 

6 

1 6 

0 

1  1 

0 

12 

9 

21 

1  96 

6 

1 

4, 

6 

1  1 

16 

1 

3 

4 

0 

31. 

6 

1 

22 

j  100. 

0 

0 

0. 

0 

16 

29. 

6 

29. 

0 

22 

1 

9 

1  6 

68. 

2 

7 

31 

8 

7 

31 

7 

1  6 

0 

32, 

6 

86. 

4 

3 

1  3 

6 

9 

33. 

9 

21. 

0 

87. 

2 

20 

90 

9 

2 

9 

1 

4 

39 

8 

19 

6 

42 

4 

18 

81  . 

8 

4 

1  8 

2 

1  6 

42. 

9 

20. 

0 

42. 

6 

16 

72 

7 

6 

27 

3 

18 

61 

6 

43 

0 

30 

3 

14 

63 

6 

8 

36 

4 

3 

64 

6 

33. 

0 

43 

9 

1  3 

69 

1 

9 

40. 

9 

1 

80. 

6  ' 

49 

0 

1 

72. 

7 

1  3 

69 

1 

9 

40. 

9 

20 

92. 

‘  i 

67  . 

1 

0 

74. 

0 

1  1 

60. 

0 

1  1 

60. 

0 

13 

116. 

0 

1  1  6 

0 

76 

0 

2 

9. 

1 

20 

90. 

9 

2 

216. 

3  1 

266 

0 

80. 

7 

3 

13 

6 

19 

86 

4 
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Figure  16  -  Actual  Target  Identification  Times,  New  Set  of  Radar  Plates 
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a  toll  of  time  and  resources  that  it  was  impossible  to  use  sufficient  plates  in 
the  validation  experiments  to  assure  a  significant  statistical  result.  The 
original  R  of  0.689  could  be  expected  to  shrink  considerably  with  new  pictures 
and  subjects  even  though  both  groups  of  subjects  had  the  same  abilities  and 
degree  of  learning.  The  original  and  validation  groups  did  have  these  quali¬ 
fications,  as  indicated  by  the  R  of  0.97  for  the  seven  common  slides;  thus  the 
obtained  shrinkage  is  attributable  to  lack  of  validity  of  the  prediction  equation 
when  utilized  with  new  imagery.  The  minimum  number  of  plates  required  to 
reject  the  null  hypothesis  for  an  R  shrunk  to  0.  Z84  would  have  been  about  110. 
This  was  determined  by  use  of  the  following  formula: 

„  ^  0. 284^  N  -  5 

4,  (N  -  5)  1  -  0. 284^  ^ 

=  0.  022(N  -  5)  . 
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LOCKDOURNE  MEAN  TARGET  -  IDENTIFICATION  TIME  (SECONDS) 


0  5  10  15  20  25  30  35  40  A5  50  55  60  65  70 

GRIFFISS  AND  PEASE  MEAN  TARGET  -  IDENTIFICATION  TIME  (SECONDS) 


Figure  17  -  Target-Identification  Times  for  Seven  Radar  Slides  Common  tfi5 
Two  Groups  of  Subjects 
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TABLE  IX 


PREDICTED  VERSUS  ACTUAL 
TARGET -IDENTIFICATION  TIMES, 
NEW  SET  OF  RADAR  PLATES 


Plate 

Mean  identification  time  (sec) 

Predicted 

Actual 

1 

12.  6 

80. 6 

2 

13.6 

216.  3 

4 

15.  7 

39.  8 

5 

8.  9 

12.  1 

6 

12.  8 

16.  0 

7 

12.  7 

31  .  7 

8 

12.  8 

7.  8 

9 

1  5.  4 

33.  9 

10 

9.  8 

10.  2 

12 

12.  4 

8.  5 

From  an  F -distribution  table,  an  F.  c.  was  sought  that  would  yield  an  F 

(In  -  I? ) 

significant  at  the  0.05  level.  For  N  =  130,  F^  ^^5  ~  2.64  from  the  above 

equation  and  F^  “  2.44  from  the  table.  For  N  =  105,  the  corresponding 

figures  were  2.20  and  2.46;  thus  the  total  plates  required  would  be  about  110. 
However,  there  appears  to  be  no  tenable  reason  to  suspect  that  with  a  large 
number  of  plates  the  correlation  would  have  been  as  high  as  0.284.  It  can  only 
be  concluded  that  no  validity  for  the  precision  equation  has  been  demonstrated. 

INFRARED  AND  OPTICAL 

Thirteen  infrared  and  15  optical  pictures  were  selected  to  check  the  application 
of  the  human-response  prediction  equation  to  these  two  types  of  imagery.  The 
infrared  pictures  were  acquired  during  daylight  at  low  altitude  in  Quick  Check 
flights  in  equipment  that  utilized  a  cooled  detector  sensitive  to  near -infrared 
energy.  The  optical  pictures  were  acquired  at  1  0,  000 -  ft  altitude  by  a  T -  1  1 
camera  with  a  6-in.  lens  at  a  scale  of  1  to  20,  000.  The  target  categories  were 
baseball  diamonds,  shopping  centers,  industrial  complexes,  drive-in  theaters, 
and  stadiums. 

Positive  glass  plates  were  prepared  and  the  same  physical  characteristics 
were  measured  in  the  same  manner  as  for  the  high-resolution  radar  pictures. 
Table  X  lists  these  values  and  the  predicted  target-recognition  times  obtained 
by  substituting  them  into  equation  Z7.  Plate  7  is  shown  in  three  parts  because 
it  contained  three  targets. 
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TABLE  X 


METRIC  I  APPLIED  TO  INFRARED  AND  OPTICAL 
PICTURES  TO  PREDICT  TARGET -IDENTIFICATION  TIMES 


Ratios  ot  transmissivity 
characteristics 

Ratios  of  differences  between  target 

Predicted 
identifica  - 

mm 

*t 

■a 

*d 

background  characteristic! 

B 

Picture 

®T 

la 

®D 

"r 

tion  time 
( sec) 

Infrared 

1 

0.  54 

1.61 

1.04 

1.03 

0.  3919 

0.2311 

0.0015 

0. 0009 

13.  1 

4 

0.90 

0.91 

1.  31 

1.  32 

0.  01  1  1 

0. 0089 

0.0738 

0.0755 

16.9 

5 

0.73 

1.  35 

1.  16 

1.  13 

0.  0999 

0.0907 

0.0221 

0.0150 

14.9 

6 

0.  34 

1.03 

0.  78 

0.90 

1.281 

0. 0009 

0.0625 

0.0111 

15.  3 

7a 

0.  86 

1.85 

1.13 

1.  14 

0.228 

0.  3905 

0.0150 

0.0172 

15.0 

7b 

0.  68 

1.02 

0.84 

0.81 

0.  1506 

0.0005 

0.0305 

0.0445 

16.4 

7c 

0.  32 

1.59 

0.  89 

0.85 

1.445 

0. 1057 

0.0136 

0.0265 

14.  3 

10 

0.  34 

1.05 

0.  78 

0.  73 

1. 281 

0. 0024 

0.0621 

0.0999 

16.  1 

1  1 

1.62 

0.65 

1.  35 

1.  31 

0.  2346 

0. 1885 

0. 1885 

0.0734 

15.7 

12 

1.03 

0.  32 

1.43 

1.42 

0.0009 

1.445 

0.2272 

0.2369 

17.  3 

13 

0.69 

0.62 

1.03 

0.98 

0.  1  393 

0. 2329 

0.0009 

0.0004 

12.9 

14 

0.  22 

2.47 

0.  88 

0.  82 

2.  765 

0.8745 

0.0164 

0.0052 

13.5 

15 

0.93 

0.98 

1.  52 

1 . 43 

0.  0053 

0.0052 

0. 1779 

0. 1293 

17.8 

16 

1. 47 

3.  52 

1.47 

1. 57 

0.  1503 

1.804 

0.  1503 

0.  2069 

15.8 

19 

1.03 

3.  47 

1. 29 

1.  36 

0.0009 

1 . 758 

0.0652 

0.953 

16.5 

Optical 

20 

1.31 

0.  52 

1.11 

1.87 

0.0734 

0.4431 

0.0109 

0.  4048 

16.0 

22 

1.40 

4.  18 

1.  36 

1.  39 

0.  1  143 

2.419 

0.0953 

0. 1094 

15.4 

23 

1.17 

1.49 

1.09 

1.08 

0.0247 

0.1611 

0.0074 

0.0059 

14.  5 

24 

1.48 

3.05 

1.53 

1. 58 

0.  1557 

1.  378 

0.  1836 

0.2129 

15.9 

25 

1. 78 

10.  65 

1. 72 

1.94 

0. 3418 

8.734 

0. 3014 

0. 4555 

15.8 

26 

1.86 

1. 84 

1.44 

1.64 

0.  3976 

0.  3835 

0.  1 344 

0. 2498 

15.6 

27 

1.18 

1.82 

1.25 

1.24 

0.0275 

0. 3695 

0.0500 

0.0465 

15.6 

28 

1.51 

2.50 

1. 37 

1.45 

0.  1723 

0.9000 

0.0999 

0.  1  397 

15.7 

29 

1.09 

0.02 

0.99 

1.05 

0.0074 

48.02 

0.0001 

0. 0238 

13.6 

30 

1. 04 

1.24 

1.01 

0.94 

0.0015 

0.0470 

0.0001 

0. 0038 

14.  3 

33 

1.15 

3.23 

1.21 

1.21 

0.0196 

1.540 

0. 3645 

0.0365 

15.2 

34 

1. 55 

3.  19 

1. 52 

1.61 

0.  1952 

1.504 

0. 1779 

0.231  1 

15.9 

35 

1.  36 

1.  30 

1.20 

1.  35 

0.0953 

0.0692 

0.0333 

0.0907 

14.9 

36 

1.  18 

1. 70 

1.  14 

1.11 

0.0275 

0.2882 

6  «  1  7  ' 

0.0109 

14.8 

37 

0.93 

2.63 

1.01 

1.01 

0.0053 

1.010 

3  -Oc 

0.0001 

13.9 
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The  validation  tests  for  the  near -infrared  and  optical  pictures  were  conducted 
with  21  interpreters  at  the  Photo  Intellij'cnce  School  at  Lowry  Air  Force  Base, 
Colo.  Glossy  prints  were  used.  The  targets  were  cut  from  duplicate  prints 
and  fixed  to  cards  with  transparent  tape.  Except  that  the  tape  darkened  them, 
the  targets  were  exact  duplicates  of  those  that  the  subjects  were  to  identify. 

A  standard  set  of  instructions  was  read  to  each  subject  (see  appendix  V).  Six 
practice  trials  were  run  and  the  tests  proceeded  in  the  same  manner  as  the 
original  tests  in  the  radar  experiments.  The  operator  timed  the  responses 
with  a  decimal  stop  watch.  These  data  are  plotted  in  figure  18  and  tabulated 
in  table  XI. 

Table  XII  gives  a  direct  comparison  of  the  predicted  and  actual  target  identi¬ 
fication  times.  The  correlation  coefficient  was  0.040  for  the  near  -  inf  ra  red 
pictures  and  0.068  for  the  optical.  By  F-test,  neither  was  significantly  dif¬ 
ferent  from  zero.  As  in  the  case  of  the  high- resolution  radar  pictures,  the 
number  of  plates  was  insufficient  to  ensure  significant  statistical  results,  but 
it  was  obvious  from  tiie  very  low  correlations  that  the  human -response  pre¬ 
diction  expression  (equation  27)  was  not  applicable  to  infrared  and  optical 
pictures. 


Figure  18  -  Actual  Target -Identification  Times,  Infrared  and  Optical  Pictures 
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TABLE  XI 


ACTUAL  TARGET -IDENTIFICATION 
TIMES,  INFRARED  AND  OPTICAL  PICTURES 


Picture 

Identification  time 
(hundredths  of  a  second) 

Number 

correct 

Percent 

correct 

Number 

incorrect 

Percent 

incorrect 

Mean 

Median 

Standard 

deviation 

Infrared 

IQ 

7.  2 

6.  0 

5.  2 

21 

100.  0 

0 

0.  0 

13 

8.  4 

6.  0 

6.  8 

21 

100.  0 

0 

0.  0 

15 

12.  4 

8.  0 

9.  6 

21 

100.  0 

0 

0.  0 

12 

13.  8 

11.0 

9.6 

20 

100.  0 

0 

0.  0 

6 

16.  3 

14.  0 

9.7 

21 

100.  0 

0 

0.  0 

10 

26.  1 

8.  5 

26.  9 

20 

95.  2 

1 

4.  8 

14 

34.9 

21. 0 

33.  2 

21 

100.  0 

0 

0.  0 

1 

36.  3 

16.  0 

63.  0 

21 

100.  0 

0 

0.  0 

11 

61. 0 

31. 0 

91  0 

21 

100.  0 

0 

0.  0 

4 

87.  1 

59.  0 

114.  6 

16 

76.  2 

5 

23.  8 

5 

119.  2 

70.  5 

103.  0 

20 

95.  2 

1 

4.  8 

16 

189.  2 

109  5 

145.  7 

10 

47.6 

11 

52.4 

7a 

9.  0 

3.  0 

20.  1 

21 

100.  0 

0 

0.  0 

7b 

11.7 

10.  0 

4.  7 

21 

100.  0 

0 

0.  0 

7c 

20.  1 

17. 0 

12.7 

21 

100.  0 

0 

0.  0 

Optical 

27 

14.9 

9.  5 

17.9 

20 

100.  0 

0 

0.  0 

24 

19.7 

13.  0 

18.  1 

20 

100.  0 

0 

0.  0 

30 

24.2 

13.  0 

30.  7 

20 

100  0 

0 

0.  0 

36 

26.4 

20.  5 

19.6 

20 

100.  0 

0 

0.  0 

25 

29.  3 

29.  5 

17.  5 

20 

100  .  0 

0 

0.  0 

33 

44.  5 

47.  0 

33.  1 

20 

100.  0 

0 

0.  0 

35 

56.  5 

44.  0 

48.  5 

20 

100.  0 

0 

0.  0 

29 

57.  3 

35.  0 

52.  8 

20 

100.  0 

0 

0.  0 

20 

59.7 

23.  0 

86.  0 

19 

95.  0 

1 

5.  0 

22 

81. 9 

32.  0 

106.  5 

17 

85.  0 

3 

15.  0 

37 

108.  5 

55.  0 

104.  0 

20 

100.  0 

0 

0.  0 

23 

1 10.  1 

54.  0 

112.  5 

15 

75.  0 

5 

25.  0 

28 

113.4 

79.  5 

106.  8 

20 

100.  0 

0 

0.  0 

34 

118.6 

69.  0 

134.4 

14 

70.  0 

6 

30.  0 

26 

217.  1 

143.  5 

151. 4 

18 

90.  0 

2 

10.  0 
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TABLE  XII 


PREDICTED  VERSUS  ACTUAL 
TARGET -IDENTIFICATION  TIMES, 
INFRARED  AND  OPTICAL  PICTURES 


Future 

Mean  identification  time  (sec) 

Predicted 

Aetna  1 

Infrared 

1 

13  1 

2  1.8 

4 

16, 9 

52.  3 

5 

14.  9 

7  1,5 

6 

15.  3 

9.  8 

7  a 

1  5.  0 

5  4 

7  b 

16.4 

7  .  0 

7  c 

14.  3 

12.  i 

10 

16.  1 

15.7 

1 1 

15.  7 

36.  6 

12 

17.3 

8  3 

1  i 

12.9 

5  0 

14 

13.5 

20.  9 

1  5 

17.8 

7.  4 

16 

15.8 

113.5 

19 

’.0.5 

4.  3 

OjJt  u  a  1 

20 

16.  0 

cc 

22 

15  4 

49.  1 

23 

14  5 

66 .  1 

24 

15,9 

11.8 

25 

15.8 

17.6 

26 

1  5  0 

1  30.  3 

27 

1  5.  6 

8.  0 

28 

15.7 

68.  0 

29 

34.  4 

13.6 

30 

14.  3 

19.5 

1  3 

15.  2 

Zb  7 

34 

15.9 

7  1.2 

35 

14  9 

33.  9 

36 

14  8 

15  8 

37 

13  9 

65.  1 
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SECTION  VI 


CONCLUSIONS  AND  RECOMMENDATIONS 


CONCLUSIONS 

Any  correlation  coefficient  obtained  as  a  measure  of  the  goodness  of  fit  of  a 
prediction  equation  to  the  target-recognition  times  used  to  establish  the  con¬ 
stants  for  the  equation  can  be  expected  to  shrink  when  applied  to  new  pictures. 
Metric  I,  which  was  developed  and  tested  in  this  study,  had  an  original  coef¬ 
ficient  of  0.689  which  shrank  to  0.284  when  applied  to  new  pictures.  This 
was  too  small  to  be  of  practical  value.  While  disappointing,  this  was  not  en¬ 
tirely  unexpected.  It  was  due  partly  to  the  unexpected  correlation  between 
the  four  variables  chosen  for  measurement.  A  contributing  factor  was  the 
difficulty  in  calculating  the  area  of  each  cell  occupied  by  a  target.  These 
calculations  had  to  be  carried  out  manually  and  in  general  were  based  on  ap¬ 
proximations. 

Had  time  permitted,  only  two  of  the  four  selected  picture  characteristics 
would  have  been  retained:  the  average  transmissivity  and  detail.  The  latter 
would  have  been  redefined  as  the  number  of  crossings  of  the  average  trans¬ 
missivity  level  by  the  signal  instead  of  as  the  absolute  derivative  of  the  signal. 
This  would  make  the  two  independent  of  each  other. 

In  addition,  at  least  one  more  characteristic  would  have  been  introduced  by 
making  use  of  the  crosscorrelation  of  the  target  with  the  whole  picture  includ¬ 
ing  the  target  as  described  in  Metric  U.  The  outcome  of  such  crosscorrelation 
is  a  surface  having  a  maximum  at  the  location  of  the  target  and  many  secondary 
peaks  at  locations  where  picture  elements  appear  more  or  less  similar  to  those 
of  the  target.  The  number  and  relative  strength  of  these  peaks  could  be  com¬ 
bined  into  one  additional  characteristic.  Also,  other  promising  data  from 
Metrics  U  and  lU  could  be  combined  with  the  foregoing  quantities  to  quantify 
image  complexity  by  automatic  devices  in  a  more  effective  manner. 

RECOMMENDATIONS 

The  development  of  a  new  metric  is  recommended.  Most,  if  not  all,  of  the 
parameters  have  been  identified  and  the  techniques  for  measuring  them  have 
been  developed.  The  slide  projector -viewer  has  been  modified  for  airborne 
use  and  currently  is  being  used  to  obtain  operator  performance  data  on  another 
program.  These  data  could  be  used  for  both  metric  derivation  and  validation 
studies. 

Some  minor  modifications  to  the  flying  spot  scanner -analog  computer  system 
and  procedures  would  increase  the  accuracy  of  the  parameter  measurements 
and  reduce  the  time  required  per  measurement.  This  would  require  adding  an 
analog-to -digital  converter  and  card-punch  unit  to  the  equipment  readout  and 
using  larger-scale  pictures  for  the  target  area  measurements. 

The  new  prediction  equation  should  be  developed  with  Metric  1  as  the  basis. 

The  average  transmissivity  variable  should  be  retained  along  with  the  revised 
detail  variable,  and  combined  with  the  most  promising  Metric  II  and  III  data 
to  form  an  equation  that  would  be  more  effective  in  quantifying  image  com- 
ple:;ity  by  automatic  devices. 


45 


appf:ndix  I 


METRIC  U,  STATISTICAL  DESCRIPTION  OF  TARGET -IDENTIFICATION 
TIME  AND  PROBABILITY  OF  CORRECT  IDENTIFICATION 

INTRODUCTION 

One  of  the  objectives  of  the  image -quantification  program  was  to  provide  a 
quantitative  analytical  description  or  model  of  the  photo -interpretation  task. 
The  photo -interpreter  was  required  to  examine  the  photograph  and  to  detect 
the  targets  of  interest  within  it.  The  prediction  of  target-identification  time 
and  the  probability  that  the  targets  would  be  correctly  detected  are  considered 
under  the  appropriate  headings  below. 

In  the  photo -interpretation  task  that  was  to  be  used  for  experimental  verifica¬ 
tion  of  Metric  II,  the  photo -interpreter  was  to  be  given,  with  each  photograph, 
a  graphic  representation  of  the  target  and  associated  cues.  The  graphic  rep¬ 
resentation  primarily  would  show  relative  locations  rather  than  shapes.  The 
examination  time  and  the  correctness  of  the  target  identification  would  be  re¬ 
corded  for  each  photograph. 

It  was  assumed  for  the  analytical  description  of  this  task  that  the  photo-inter¬ 
preter  first  would  select  a  major  reliable  cue,  and  then  coarsely  (in  accord¬ 
ance  with  the  size  of  the  cue)  scan  the  photograph  for  the  presence  of  this  cue. 
At  each  scan  point  that  appeared  to  contain  the  cue,  the  photograph  would  be 
examined  in  greater  detail  with  a  fine  scan  for  the  detection  of  the  actual  tar¬ 
get  or  other  cues  that  would  reinforce  the  reliability  of  the  target  location. 

The  probabilistic  determination  of  the  coarse  and  fine  scanning  times  for  a 
given  photograph  would  provide  an  estimate  of  the  time  required  for  proper 
examination  of  the  photographs.  These  times  were  taken  to  be  functions  of 
(1)  target  complexity,  (2)  the  similarity  between  the  targets  and  their  immedi¬ 
ate  environments,  (3)  the  number  of  false  match  points,  and  (4)  the  number  of 
cues  required  to  identify  the  target  and  their  cue -to -environment  similarity. 
The  target  parameters  were  to  be  evaluated  in  terms  of  statistical  data  ob¬ 
tained  through  television  raster  scans  of  the  photographs  restricted  to  targets 
having  the  same  brightness  as  the  cues.  The  video  signal  output  of  the  scanner 
was  to  be  fed  through  a  low-pass  filter  to  limit  the  effective  resolution  of  the 
signal  to  approximately  one -fifth  of  the  cue  length.  Analytical  relationships 
between  these  experimentally  obtained  distribution  functions  and  the  scanning 
times  would  determine  the  desired  prediction  model  of  the  photo -interpretation 
time . 

A  photo -interpreter  recognizes  a  target  primarily  by  its  shape  or  pattern,  as 
in  airfield  runways,  railroad  yards,  and  petroleum  tank  farms,  or  by  its  tex¬ 
ture  and  intensity,  as  in  cities  and  bodies  of  water.  These  target  character¬ 
istics,  in  addition  to  the  relative  locations  of  targets,  may  be  presented  to  the 
photo -interpreter  as  part  of  his  briefing;  more  usually  they  are  acquired  by 
him  during  previous  photo -interpretation  training  and  experience.  Therefore, 
the  probability  that  a  photo -interpreter  will  correctly  identify  the  targets  of  a 
given  picture  is  essentially  the  probability  that^these  unique  target  character¬ 
istics  are  indeed  associated  with  the  targets  of  the  picture  and  are  absent  from 
the  background. 
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The  evaluation  of  this  probability  of  the  presence  of  target  recognition  cl^arac- 
teristics  requires  that  these  characteristics  be  converted  to  a  set  of  target¬ 
defining  parameters  that  can  be  experimentally  obtained  from  the  picture.  In 
practice,  the  picture  is  subdivided  into  a  matrix  of  equal  target-size  areas, 
and  a  set  of  parameters  is  obtained  for  each  target-size  area.  Each  target 
class  that  is  of  interest  to  the  photo-interpreter  can  be  described  in  terms  of 
a  unique  set  of  experimentally  acquired  statistical  density  functions  for  the 
distributions  of  these  parameters.  Insertion  of  the  parameter  values,  meas¬ 
ured  for  a  given  target-size  area,  into  the  target-class  statistical  density 
functions  permits  evaluation  of  the  probability  that  this  set  of  parameters, 
and  therefore  the  given  area  that  the  set  represents,  is  associated  with  a 
certain  target-class. 

TARGET -IDENTIFICATION  TIME 

Qualitative  Description  of  Photo -Interpretation  Task 

The  photo -interpreter  would  be  given  a  graphic  representation  (primarily  with 
respect  to  location  rather  than  shape)  of  each  target  and  its  associated  cues. 
With  this  representation  he  would  be  required  to  examine  the  photograph  and 
correctly  detect  the  target. 

Many  targets,  such  as  airfields,  industrial  sites,  military  installations,  and 
railroa*..  marshalling  yards,  could  not  be  used  as  primary  cues  because  they 
were  too  small  for  the  photographic  scale  (20  by  20  nautical  miles).  Therefore, 
large  neighboring  cues,  such  as  culture  and  rivers,  and  smaller  reinforcing  cues, 
such  as  roads,  railroads,  and  drainage  patterns,  were  to  be  used  by  the  pho¬ 
to-interpreter  to  pinpoint  the  locations  of  the  targets.  In  the  photo -interpreta - 
tion  procedure,  a  large  and  reliable  cue  would  be  selected  and  a  coarse  scan 
(the  coarseness  of  which  would  be  set  by  the  size  of  the  cue)  would  be  used  to 
locate  and  identify  the  selected  cue.  To  obtain  maximum  simplicity  and  re¬ 
liability  in  locating  the  target,  the  photo-interpreter  would  select  his  primary 
cue  to  be  the  brightest  (such  as  cities  or  culture)  or  the  darkest  (such  as  rivers 
or  lakes)  pattern  in  the  photograph.  The  photo -interpreter  then  would  be  con¬ 
cerned  only  with  those  patterns  having  the  same  intensity  level  as  the  given 
cue. 

Once  the  cue  had  been  located,  a  more  detailed  search  in  the  neighborhood  of 
the  cue  would  be  necessary  to  verify,  through  the  remaining  auxiliary  cues 
given  in  the  reference  graphic,  that  the  located  cue  was  indeed  the  true  one. 

The  complete  interpretation  process  therefore  would  be  one  of  (1)  scanning 
for  patterns  of  the  same  brightness  level  as  the  cue,  (2)  locating  patterns  that 
matched  the  shape  of  the  cue,  and  (3)  determining  that  the  cue  was  either  true 
or  false.  Statistical  estimates  of  the  times  required  for  each  of  the  three  sub¬ 
processes  are  given  in  the  following  paragraphs.  The  terminology  used  to  ob¬ 
tain  the  statistical  time  estimates  is  given  in  table  XUJ. 

Location  of  Targets 

The  initial  ccarse  scan  to  locate  the  targets  of  proper  brightness  (referred  to 
herein  as  cue -level  or  i-level  patterns)  will  become  more  difficult  and  time- 
consuming  as  the  similarity  between  the  intensity  levels  of  the  target  and  its 
immediate  en/ironment  increases.  A  measure  of  patte rn -to -environment  dis¬ 
similarity  may  be  obtained  from  the  statistical  distribution  of  the  values  for 
the  gradients  that  exist  at  the  ends  of  the  i-lcvel  patterns. 
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TABLE  Xm 


TERMINOLOGY  FOR  TIME  FOR 
CORRECT  TARGET  IDENTIFICATION,  METRIC  II 


C  Photo-interpreter  capacity  (bits  per  second)  associated  with  his  rocog- 

^  nition  of  the  cue 

Capacity  of  photo-interpreter  (bits  per  second)  to  scan  a  sequence  of 
nonuniform  targets 

Number  of  different  i-level  target  lengths  in  a  given  picture 

Absolute-value  intensity  gradient  before  and  after  each  i-level  target 
Ij  Cue  complexity  (bits)  represented  by  information  transmitted  by  cue 

Number  of  cues  associated  with  the  desired  target 

Nj  Number  of  false  match  points  obtained  up>on  the  use  of  a  given  cue 

N^  Total  number  of  subareas  examined  in  a  given  picture 

N(x^)  Number  of  i-level  targets  with  a  length  of  cells 

Probability  of  encountering  an  i-level  cell 

P^  Probability  that  the  matched  cue  is  not  associated  with  the  desired 
target  (probability  of  a  false  match  point) 

P^  Probability  that  a  cue  is  matched 

T  Expected  total  time  required  by  photo-interpreter  to  locate  the  de¬ 

sired  targe*  within  the  photograph 

T^  Total  time  used  by  photo-interpreter  to  locate  a  target 

T^  Total  time  required  to  distinguish  between  true  and  false  target  areas 
when  all  the  F>attern  cues  are  used 

tg(g)  Average  time  delay  associated  with  extraction  by  the  photo-interpreter 
of  a  target  with  edge-gradient  g 

T  Average  photo-interpreter  time  delay  resulting  from  the  target-to- 

environment  similarity 

x^  Lengths  of  i-level  targets 

At  Time  required  to  distinguish  between  the  true  and  false  target  areas 
when  a  single  cue  is  used 
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Experimentally,  the  instantaneous  absolute -value  gradient,  g,  will  be  avail¬ 
able  before  and  after  each  i -level  pattern;  however,  the  average  gradient  at 
each  edge  will  be  a  better  measure  of  target -to -environment  similarity  and 
should  be  used  when  possible.  A  probability  density  function  of  the  gradient, 
g,  can  be  derived  from  its  experimentally  obtained  first  two  mo  nents  (ref¬ 
erence  7).  It  will  be  a  truncated  Gaussian  (reference  8): 

P(g)  =  exp  -  (a.  +  b.g  +  c.g^)  ,  (28) 


where  g  >  0. 


The  term  t  (g)  represents  the  time  delay  associated  with  the  locating  of  a  tar¬ 
get  with  edge -gradient  g.  The  average  (over  the  entire  scene)  time -delay  per 

target,  T  ,  because  of  target-to-environmental  similarity,  will  become 
8 


1' 

s 


/"  l-)[ 
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exp  -  (a.  +  b.g  +  c.g  ) 


dg 


(29) 


A  secondary  time  delay  in  scanning  a  sequence  of  i -level  targets  (those  having 
the  same  brightness  as  the  major  cue)  will  be  incurred  because  of  nonuniform¬ 
ity  of  target  lengths.  Such  nonuniformity  can  be  measured  from  the  informa¬ 
tional  entropy  of  the  distribution,  p(xp,  of  i -level  target  lengths: 

X.  =  00 
1 

P(x.)  log2  P(x.)  . 

X.  =  1 
1 


in  bits  per  target.  If  represents  the  capacity  of  the  photo -interpreter  in 
bits  per  second  to  scan  a  sequence  of  nonuniform  targets,  the  average  time 
delay  per  target  because  of  target  nonuniformity  will  be 


x.(max) 


X.  =  1 


1 


P(x.)  log^  P(>:-)  • 


(30) 


The  total  number  of  i -level  targets  in  the  photograph  will  be 


nP. 


<x.> 


where 
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n  =  the  number  of  cells  in  the  photograph  (cell  length  is 
assumed  to  be  one -fifth  the  length  of  cue), 

P.  =  probability  of  encountering  an  i -level  cell,  and 
<x^>  =  average  length  of  i -level  targets. 


The  total  target  location  time,  T/.  will  be 
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P(x.)  log^  P(x.)  . 


Scanning  of  True  and  False  Match  Points 


(31) 


Once  a  possible  match  is  found  for  the  primary  cue  of  the  reference  graphic, 

the  photo -interpreter  will  initiate  a  finer  and  more  detailed  scan  to  locate  and 

identify  the  remaining  cues  of  the  graphs  and  thus  verify  that  the  possible 

match  point  is  either  true  or  false.  The  number  of  false  match  points  that 

probably  will  be  encountered  in  a  given  scan  will  be  obtained  by  assuming  that 

the  cue  is  composed  of  a  group  of  one -dimensional  patterns  having  C  different 

lengths  and  that  there  are  N(x.)  targets  within  a  given  x.  length  class.  The 

probability,  P  ,  then  will  be 
o 


p 

o 


N(x  ) 

P(x.) 


(32) 


P(x.),  the  distribution  in  the  photograph  of  the  i-level  target  lengths,  is  pre¬ 
sented  in  connection  with  aquation  31. 


The  deviation  from  perfect  match  will  be  assumed  to  have  taken  place  inde¬ 
pendently  at  any  of  the  cells  of  the  cue  so  that,  for  targets  that  are  almost 
matched,  the  probability  of  cue  match,  P  ,  can  be  expressed  as  a  binomial 
distribution  in  terms  of  Pa^,  the  probability  that  an  individual  cell  of  the  cue 
is  matched.  The  cells  will  be  assumed  to  have  lengths  equal  to  one -fifth  that 
of  the  cue.  If  there  are  N  cells  in  the  cue,  the  value  for  P  will  be  obtained 
from  ^ 


N 


P  = 
o 


i  =  C  N(x.) 

=  n  P(x.)  '  .  (33) 

i  =  1 
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Since  the  briefing  knowledge  concerning  the  primary  cue  will  not  necessarily 
be  complete,  the  actual  matching  of  the  primary  cue  with  its  counterpart  in 
the  scene  will  not  be  perfect.  The  degree  of  mismatch  will  depend  upon  the 
cue.  Bodies  of  water,  since  they  are  usually  adequately  surveyed  and  there¬ 
fore  accurately  located,  will  exhibit  a  high  degree  of  match;  culture,  con¬ 
versely,  cannot  be  accurately  predicted  and  will  in  general  show  a  low  de¬ 
gree  of  match.  If  there  are  N  cells  in  the  cue  and  if  m/N  represents  the 
degree  of  match  between  the  primary  cue  and  its  counterpart  in  the  scene,  the 
probability  of  match,  P  ,  at  other  points  in  the  scene  (which  in  fact  is  the 
probability  that  a  false  match  point  has  been  encountered)  will  be 
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(34) 


The  problem  of  establishing  a  correspondence  or  match  between  each  cue  may 
be  thought  of  as  matching  a  square  area  just  large  enough  to  include  the  cue. 
The  quantitized  possible  match  positions  are  referred  to  as  a  correlation 
matrix. 

If  N  is  the  number  of  points  in  this  correlation  matrix  and  L  (in  cells)  is  the 
average  length  of  the  cue,  the  probable  number  of  false  match-point  loops  en¬ 
countered  in  the  scan  will  be  approximately 

N  P  P. 

* .  _  m  m  1 

Nf  =  - E - ■ 


The  time,  At  ,  spent  in  the  investigation  of  a  false  or  true  match  point  will 
be  a  function  oV  the  number  and  complexity  of  the  reference  cues  that  must  be 
identified  before  the  decision  is  made  that  the  match  point  is  true  or  false. 

A  cue -identification  time,  t.,  may  be  assigned  to  a  given  class  of  cues  so  that 
the  total  match-point  time  (Kt  )  becomes 
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j  =  1 
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(36) 


The  value  for  j  will,  of  course,  depend  upon  the  number  of  cues  required  to 
render  the  true  or  false  match-point  decision.  In  general,  all  the  cues  pres¬ 
ent,  j  =  jj,  in  the  reference  graphic  will  be  used  in  the  true  match-point 
decis^n,  while  some  lesser  jf  (probably  2  or  3)  will  provide  a  reliable  false 
match-point  decision.  With  equation  34  used  to  evaluate  P  ,  the  total  scan 
time,  T^,  for  the  analysis  of  the  true  and  false  match  poi^s  will  become 
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J, 


t. 

J 


(37) 


Cue  Examination  Time 


The  tim  required  to  examine  a  cue  depends  upon  the  complexity  of  its  pattern. 
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Cue  patterns  vary  from  the  very  simple  straight  roads,  railroads,  or  circular 
blobs  that  represent  small  cultural  areas  to  the  more  complex  drainage  patterns, 
meandering  rivers,  and  large  cultural  areas. 

The  complexity  of  a  cue  may  be  quantitatively  determined  by  calculating  the 
amount  of  information,  I,  transmitted  by  the  cue  (reference  9)-  f'or  example, 
the  amount  of  information  transmitted  by  an  auxiliary  "straight-road"  cue  due 
to  its  orientation  and  location  may  be  calculated  as  follows.  The  probability  of 
a  certain  orientation  of  the  road  is  l/36.  This  assumes  that,  the  photo-inter¬ 
preter  can  discern  the  orientation  of  the  road  within  a  10-deg  azimuthal  inter¬ 
val.  The  probability  that  the  cue  is  in  a  certain  location  will  be  the  ratio  of  the 
auxiliary  cue  length  (for  the  road  this  is  approximately  200  ft)  to  the  length  of 
the  primary  cue  (assumed  to  be  10,  000  ft),  and  for  the  road  will  be  equal  to 
l/50.  The  information  transmitted  by  the  road  will  be 

Kroad)  =  -  log^  ^  -  log^  5^  , 


=  9.07  bits  . 


(38) 


In  this  manner,  one  can  obtain  a  quantitative  estimate,  in  bits,  of  the  complex¬ 
ity  associated  with  the  various  cues  of  a  given  graphic.  The  actual  measured 
times  required  by  the  photo-interpreter  to  recognize  the  various  cues  will  pro¬ 
vide,  in  conjunction  with  the  above  "informational"  values  of  cue  complexity, 
an  estimate  of  the  photo-interpreter  capacity,  Cp,  in  bits  per  second  to  recog¬ 
nize  cues.  The  average  human  information  rate  is  approximately  25  bits  per 
second  (references  10  and  ll)and  thus  it  appears  that  a  cue  examination  time  of 
leas  than  a  second  will  be  sufficient.  However,  the  actual  recognition  time 
per  cue  probably  will  be  much  longer  because  the  photo -interpreter  will  usu¬ 
ally  refer  back  and  forth  from  the  photograph  to  the  reference  graphic  during 
the  fine  scanning  process  for  reinforcement  of  his  stored  description  of  the 
cue.  Furthermore,  a  discriminating  task,  wherein  high  accuracy  rather  than 
speed  is  the  goal,  is  almost  always  more  difficult  (reference  10,  pp  130  and 
131). 


The  above  calculations  of  cue -recognition  time  assume  that  the-  contrast  be¬ 
tween  the  brightness  of  the  cue  and  its  immediately  adjacent  background  (the 
edge  gradient)  will  be  sufficiently  high  so  that  no  time  delay  will  be  incurred 
in  the  reading  of  the  cue  pattern.  If  the  contrast  is  low,  an  additional  time- 
delay  term  will  be  required.  This  term  is  of  the  order  of  t  (g).  Thus,  the 
cue -examination  time,  tj,  for  a  given  j-cue  will  be  a  function  of  the  complex¬ 
ity  of  the  j-cue,  I.;  the  photo -interpreter  capacity,  Cp',  and  if  applicable,  the 

time  delay,  t  (g.)"!  associated  with  the  contrast,  g.,  of  the  given  cue,  or 
s  J  J 


Total  Photo  Examination  Time 


(39) 


The  total  time,  Tg,  for  examination  of  the  photograph  will  be  the  sum  of  the 
expected  pattern -location  and  match  point  analysis  times  (equations  31  and 
37): 
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Equation  40  will  give  an  average  value.  Since  the  total  time  will  be  made  up 
of  several  independent  contributions,  the  variance  will  be  the  sum  of  the  in¬ 
dividual  variances  of  the  several  contributions.  The  various  terms  of  the  ex¬ 
pression  will  be  either  picture -derived  and  require  statistical  readout  of  the 
photograph  and  reference  graphic  or  will  depend  upon  the  photo -interpreter 
and  be  obtained  from  actual  photo -interpretation  tests. 


Measurement  of  Picture  Parameters 


Experimental  Acquisition  of  Picture-Derived  Terms 


The  picture -derived  terms  will  apply  to  those  patterns  that  have  the  same 
brightness  as  the  primary  cue  and  will  include  the  two  distribution  functions, 
P(g)  and  P(x.),  the  pattern  probability,  and  average  length,  P.  and  x..  The 
average  cue  'length,  L,  the  number  of  cells  in  the  picture,  n,  'the  number  in 
the  cue,  N,..,  and  the  number  in  the  equivalent  correlation  matrix,  N  ,  will 
be  obtained  from  the  picture.  Also  required  will  be  an  automatic  or  visual 
observation  of  the  number  and  complexity  of  the  reference  cues,  and  a  re¬ 
stricted  scan  of  the  target  area  on  the  problem  scene  to  obtain  the  average 
cue -to -background  contrasts.  These  measurements  will  permit  calculation 
of  the  cue -recognition  times,  t ^  . 

The  majority  of  the  photographic  data  will  be  derived  from  two  quadrature 
raster  scans  of  the  photograph.  For  the  initial  coarse  scan,  the  photo-inter¬ 
preter  will  adjust  the  coarseness  of  the  scan  and  therefore  the  cellular  reso¬ 
lution  of  his  field  of  vision  to  accommodate  the  size  of  the  cue.  The  cellular 
size  of  the  two  quadrature  data -acquiring  scans  should  be  made  equal  to  that 
used  by  the  photo -interpreter  in  his  coarse  scan  of  the  photograph  (approxi¬ 
mately  one-fifth  the  average  cue  length).  In  most  cases,  this  will  require  that 
the  video  from  the  flying  spot  sc^lnner  pass  through  a  low-pass  filter  with  the 
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filter  time  constant  being  determined  by  the  cue  size.  In  addition,  a  power - 
relationship  with  an  exponent  typically  equal  to  one-half  will  exist  between  the 
physical  brightness  of  a  given  image  and  psychophysical  judgments  of  that 
brightness  (reference  12).  The  video  from  the  flying  spot  scanner  should  be 
further  modified  to  provide  this  same  one-half  power  relationship.  All  data 
measurements  are  to  be  ma.de  with  this  modified  video. 

To  simplify  the  target  extraction  process  and  make  it  more  reliable,  the  pho¬ 
to-interpreter  will  select  fDr  his  primary  clue  a  prominent  target  that  is  either 
the  brightest  or  the  darkest  in  the  picture.  If  eight  distinguishable  brightness 
bands  characterize  a  given  photograph,  the  picture  patterns  of  interest  will  lie 
in  either  the  top  or  the  bottom  brightness  band.  However,  a  photo -interpreter 
can  distinguish  a  much  finer  brightness  difference  than  that  accorded  by  eight 
bands.  This  will  occur  when  two  narrowly  separated  levels  of  brightness  as¬ 
sociated  with  a  fine  brightness -difference  characterize  adjacent  patterns.  Ac¬ 
cordingly,  a  pattern  lengtli  x.  will  be  of  interest  when  the  pattern  brightness 
lies  in  the  same  brightness  band  as  the  cue  and  the  change  in  brightness  be¬ 
tween  the  target  and  its  neighbor  at  either  edge  of  the  length  exceeds  an  em¬ 
pirically  obtained  minimum  distinguishable  value,  g  .  . 

min 

The  measurement  and  recording  of  the  above  pattern  lengths  will  permit  deri¬ 
vation  of  an  analytical  distribution  function,  P(x.),  of  these  lengths.  Integra¬ 
tion  of  all  the  rrieasured  lengths,  x^,  will  permit  the  evaluation  of  Pj  (the  proba¬ 
bility  that  an  xj  pattern  exists).  The  distribution  function,  P(g),  may  be  de¬ 
rived  in  a  similar  manner  through  the  measurement  and  recording  of  the 
gradients  that  occur  at  the  edges  of  the  pattern  lengths,  Xj.  The  remaining 

parameters  n,  N  ,  and  N  will  be  a  function  of  the  cell  size,  which  is  assumed 
*  c  m 

to  be  one -fifth  of  L,  the  average  length  of  the  cue. 

Empirical  Photo -Interpretation  Terms 

The  remaining  terms  such  as  the  pattern  location  time,  tg(g),  the  capacity  of 
the  photo-interpreter  to  discern  nonuniform  patterns,  Cj^,  i  id,  in  a  partial 
sense,  the  cue  identification  times,  tj,  will  be  obtained  from  actual  photo - 
interpretation  tests. 

Summary 

The  above  photo -inter pretation  model  will  be  a  quantitative  estimate  of  the 
times  required  for  the  photo -interpreter  to  complete  his  several  subtasks, 
which  begin  with  a  reference  graphic  that  contains  a  complex  of  target  cues 
and  culminate  in  the  recognition  of  the  target  of  interest.  The  photo-inter¬ 
pretation  procedure  is  assumed  to  be  that  of  first  scanning  the  photograph  in 
a  coarse  manner  for  the  location  of  the  most  prominent  and  reliable  graphic 
cue.  Then  those  areas  that  are  thought  most  likely  to  contain  the  target  will 
be  examined  in  detail  through  the  use  of  the  other  graphic  cues  for  verifica¬ 
tion  of  this  target  area  and  recognition  of  the  desired  target. 

It  is  presumed  that  the  photo -interpreter  can  largely  ignore  all  targets  except 
those  that  are  at  the  same  brightness  level  as  the  target  cue.  Accordingly, 
the  model  will  use  such  experimentally  obtained  functions  of  the  cue -level 
targets  as  the  distributions  of  their  pattern  lengths  and  of  their  pattern-to- 
environment  contrasts  along  with  the  percent  of  the  photographic  area  that 
these  patterns  occupy. 
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These  will  be  combined  to  form  an  estimate  of  the  coarse  scan  time  and  of  the 
number  of  false  cue -location  points  that  are  to  be  examined  further.  The  com¬ 
plexity  of  the  environment  in  which  the  cue  is  located  as  well  as  the  complexity 
of  the  cue  itself  will  be  incorporated  into  the  estimate  of  the  fine -scan  time  re¬ 
quired  to  examine  the  smaller  areas  within  which  the  cues  and  targets  are  as¬ 
sumed  to  lie. 

PROBABILITY  OF  CORRECT  TARGET  IDENTIFICATION  BY  THE  PHOTO¬ 
INTERPRETER 

Target  Signatures 

Table  XIV  gives  the  terminology  used  in  the  derivation  of  the  probability  *^hat  the 
photo -interpreter  will  correctly  identify  the  target  area.  A  target  usually  is 
recognized  by  a  two-dimensional  uniqueness  of  either  its  pattern,  texture,  or 
intensity,  or  all  three.  The  two-dimensional  uniqueness  of  the  target  that  is 
observed  by  the  photo -interpreter  is  converted  via  a  TV  raster-like  scan  into 
a  set  of  characteristic  one -dimensional  target  signatures.  For  example,  a 
TV  scan,  when  applied  to  suitable  differentiation  (to  mark  the  pattern  edges) 
and  interval -measuring  (to  obtain  the  pattern  widths)  circuitry,  will  generate 
a  sequence  of  edge-to-edge  spacings  that  is,  in  effect,  a  characteristic  signa¬ 
ture  of  the  geometrical  uniqueness  of  the  target  pattern.  In  addition,  the  se¬ 
quential  readout  of  the  lengths  of  the  edges  forms  a  useful  pattern-uniqueness 
signature.  Restrictive  storage  of  only  those  edges  that  follow  another  by  a 
certain  allowable  distance  permits  generation  of  other  signatures  that  would 
further  enhance  the  targets'  pattern  uniqueness.  This  latter  edge  restriction 
would  permit  the  preferential  examination  cf  targets  of  a  certain  width. 

The  above  signatures  are  representative  of  the  target  only  for  a  certain  pre¬ 
scribed  target  orientation  and,  therefore,  would  be  valid  only  if  the  target 
orientation  is  known  a  priori.  In  the  general  case,  target  orientation  is  not 
known,  so  that  some  preferred  orientation  must  be  used  with  each  signature. 
Such  an  orientation  for  the  edge -length  signature  is  that  in  which  the  average 
edge  length  is  maximum;  the  preferred  edge -spacing  signature  has  a  minimum 
average  spacing.  In  addition,  the  complete  ensemble  of  signatures  that  are 
obtained  at  the  many  different  possible  pattern  orientations  forms  an  orienta¬ 
tion-sensitive  sequence  of  measurements  that  provides  useful  target -classifi - 
cation  information. 

Many  targets,  such  as  cities  and  bodies  of  water,  are  recognized  primarily  by 
their  texture  and  intensity  rather  than  their  over -all  shape.  Signatures  that 
characterize  target  intensity  and  texture  are  available  through  the  readout  of 
the  amplitude  and  derivatives  of  the  original  flying  spot  scanner  video  before 
it  is  processed  to  form  the  target  edges. 

Target  Signature  Parameters 

To  facilitate  the  experimental  processing  of  the  data  and  to  simplify  the  com¬ 
putation  of  the  required  probabilities,  the  many  measurements  of  a  target  sig¬ 
nature  sequence  will  be  reduced  to  a  much  smaller,  but  equivalent,  set  of  pa¬ 
rameters.  Each  generated  target  signature  may  be  likened  to  a  sampled 
electrical  signal  and,  as  such,  can  be  described  in  terms  of  its  energy  and 
power,  its  d-c  and  a-c  contents,  and  its  spectral  density.  If  the  attribute  to 
be  measured  at  a  given  point,  x^,  in  a  scan  is  M(x.)  and  if  during  the  course 
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Symbol 


*hi'  ^lo 


TABLE  XIV 

TERMINOLOGY  FOR  PROBABILITY  OF  CORRECT 
TARGET  IDENTIFICATION,  METRIC  U 


j  D-c  content  or  average  amplitude  of  the  high*  and  low-fre<juency  signatures,  f^^.  and 
!  Experimentally  obtained  coeffic  icnts  that  determine  the  Gaussian  probability  density  function. 

I 

'  Target  class  j 

I 

I  Average  deviation  between  the  recognition  vector  difference,  a  obtained  by  the  photo* 

I  interpreter  and  the  model -acquired  difference. 

!  Total  energy  of  f(x-) 

I  I'nergy  of  the  high*  and  low*frequenc  y  signatures,  f^^^  and  f|^ 
j  D-c  content  of  ffx.l 


fh^fx^l,  f|^(x^)  ;  High-  and  lou, -frequenc y  target  signatures  derived  from  f(x.l 


CAk).  G.fyl* 


I  The  target  signature  or  sequence  of  N  parameter  measurements 
Absolute-value  gradient  at  i*level  target  edges  (both  first  and  second  moments) 


P<A.».  P(Uj» 


P(U.|A) 


<‘*  ><v 


I  Parameter  measurement  obtained  at  the  given  sc  anning  point 

! 

'  Number  of  scanning  points  within  a  given-sise  target  area 
j  The  number  of  i-levrl-intensity  target  edges 

i  Probability  that  the  photo-interpreter  has  made  the  correct  "no-target**  decision  at  each  of  the 
j  "no-target"  subareas,  that  is.  he  has  correctly  located  the  target 

I  Experimentally  obtained  prior  probabilities  of  obtaining  target  class  A  and  cumulative  parame- 
I  ter  value  respectively  " 

U  )  I  Conditional  probability  that  a  given  subarea  represents  target  class  A.  given  the  set  of  n 
"  ^  cumulative  parameter  values  through  ^ 

Gaussian  probability  density  function  of  the  distribution  of  cumulative  parameters.  U^,  in 
*  target  c  lass  A^ 

Target  class  A^  recognition  vector  obtained  from  a  no-target  subarea  i 

'  Target  recognition  vector  for  class  A^  obtained  by  evaluating  Iog^|p(A^|U| 

The  value  of  the  target  class  A  recognition  vector  a*iSoc  lated  with  a  subarea  known  to  contain 
the  target  class  ^ 

Average  transmissivity  of  X  and  Y  scans  of  a  gi*  en  subarea 
Transmissivity  second  moments 

Experimentally  obtained  value  for  a  given  cumulative  measurement  as  f^^.  E.  ^fo*  ^hi 

E|^  of  a  giv«  *1  parameter  signature,  f(x^) 

Location  of  a  given  scanning  point 


Lengths  associated  with  i-intensity-level  targets 

Equivalent  difference  between  the  target  and  no-target  subarea  recognition  vectors  as  obtained  by 
the  photo-interpreter 


Variance  of  the  deviations  between  the  photo- interprecer-obtajned  recognition  vector  differences 
and  the  model-acqturrd  recognition  differences  (both  C  and#^  reflect  the  abilities  of  the  photo- 
interpreter  and  are  obtained  from  the  results  of  many  photo- interpretation  experiments) 


—— 1>  Average  of  absolute-value  gradients 


•  tf  hv\ 

<  *->.<“> 
dy 


Gradient  second  moments 


Cumulative  parameters  suggested  for  measurement 
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of  a  single  scan  there  are  N-  such  measurements,  the  resulting  signature  of 
the  measured  target  will  be  a  sampled  signal: 


f(x) 


=  N 


E 

i  =  1 


M(x.)  5 (x  -  X. )  , 
1  1 


(41) 


wherein  the  sampling  interval  depends  upon  the  attribute  to  be  measured  and 
may  or  may  not  be  independent  of  time.  The  d-c  content,  f^^,  and  energy,  E, 
of  the  complete  signature  will  be  obtained  by 


and 


(42) 


E  = 


(43) 


A  third  parameter  N,  the  number  of  sample  points  in  the  signature,  can  be 
obtained  readily. 

An  analysis  of  the  information  provided  by  a  parameter  in  a  recognition  decision 
indicates  that  a  computed  parameter  adds  no  more  information  to  the  decision 
than  is  already  available  from  the  individual  component  parameters  that  entered 
into  its  computation  (see  appendix  IV). 

The  parameters  E,  and  N  describe  the  signature.  Additional  energy  para¬ 

meters  that  may  be  used  to  describe  an  electrical  signal,  such  as  power,  and 
d-c  and  a-c  energies  and  amplitudes,  will  be  computed  from  the  original  three 
parameters.  Other  parameters  v/ill  be  necessary  to  describe  the  signature 
spectrally.  Such  information  is  a  function  of  the  interdependence  of  adjacent 
signature  measurements. 

Since  the  sampling  interval  is  arbitrary  and  completely  independent  of  time, 
conventional  high-  and  low-pass  filters  for  temporal  analysis  of  the  spectral 
content  of  electrical  signals  are  useless.  However,  high-  and  low-pass  filters 
are,  in  effect,  differentiating  and  integrating  circuits,  and  the  operations  of 
differentiation  and  integration  can  certainly  be  applied  to  tht.  pattern  signatures. 

Accordingly,  the  output  of  an  equivalent  high-pass  filter  will  become 

i  =  N 

fj^.  =  ^  ]  |m(x.)  -  M(x.  +  l)j|£(x  -  x.)|  ,  (44) 


and  that  of  a  low-pass  signal  will  become 
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i  -  N 

f/o  =  ^  [m(x.)  +  Mix.  +  1)]  [8(x  -  X.)]  .  (45) 

i  =  1 

The  d-c  content,  or  amplitude,  A,  and  energy,  E,  of  these  spectral  components 
(the  number  of  sample  points  N  has  already  been  explained)  will  be 

i  =  N  - 

i  =  1 
i  =  N  - 

^hi  =  E 

i  =  1 
i  =  N  - 

=  E 

i  =  1 
i  =  N  - 

=  E 

i  =  1 

Further  differentiation  and  integration  of  the  high-  and  low-pass  spectral  com¬ 
ponents  (  equations  44  and  45)  will  provide  more  detailed  spectral  descriptions 
of  the  signature.  Greater  spectral  detail  will  be  required  as  the  interdependence 
.among  adjacent  measurements  extends  beyond  one  or  two  intervals. 

Statistical  Distributions  of  Target  Parameters  in  Reference  Target  Classes 

The  previously  described  parameters,  generated  from  the  target  signatures, 
will  form  from  each  signature  a  set  of  numbers.  Up  •  •  U^^,  for  calcula¬ 

tion  of  the  desired  target  detection  probability.  However,  before  this  can  be 
done,  conditional  probability  distributions,  P(UjAj),  of  the  various  parameters 
Uj  through  Uj.j  must  be  determined  for  each  target  clase  Aj.  These  will  be  ob¬ 
tained  experimentally  as  follows.  The  prior  readout  of  many  sample  patterns 
of  a  given  target  class  will  provide  a  set  of  values  for  each  parameter.  From 
each  set  of  values,  an  experimental  probability-density  function  of  the  following 
form  for  a  single  parameter  will  be  derived: 


1 

[m(x.)  -  M(x.  +  1)]  ,  (46) 

1  2 

[m(x.)  -  M(x.  +  1)]  ,  (47) 

[m(x.)  +  M(x.  +  1)]  ,  (4  3) 

2 

[m(x.)  +  M(x.  +  1)]  .  (49) 


P(U 


A^) 


=  e 


-(A..  +  B..U.  +  C..U.  ) 


(50) 


Such  a  function,  a  truncated  Gaussian,  can  be  obtained  (see  appendix  IV)  for 
any  set  of  values  that  are  clustered  in  a  single  group,  whether  the  group  be 
tight  or  loose.  However,  value  sets  that  contain  two  or  more  distinct  clusters 
cannot  be  adequately  represented  by  the  above  probability  density  function. 
Such  data  are  to  be  divided  into  as  many  portions  as  there  are  clusters,  with 


59 


a  probability  density  function  of  the  above  form  assigned  to  each  cluster.  De¬ 
tails  of  the  conversion  of  both  single  and  multicluster  data  sets  to  the  above 
exponential  probability  density  functions  appear  in  appendix  IV. 


Target  Detection  Probability 


Target  Recognition  Vector 


In  the  evaluation  of  the  target  recognition  vector  from  which  the  target  detection 
probability  will  be  obtained,  use  will  be  made  of  the  conditional  probability, 
P(Aj|Uj  .  .  .  U^),  that  a  given  target-size  area  tnat  is  characterized  by  a 
parameter  set  Uj  .  .  .  belongs  to  a  target  class  Aj.  In  accordance  with 
Bayes'  rxile  (reference  13,  p  318),  and  since  the  parameters  will  be  inde¬ 
pendent  of  one  another,  the  probability  of  target  class  Aj  can  be  placed  in 
terms  of  the  experimentally  obtained  parameter  density  function  P(Ui|Aj),  or 

P(A.|U,  .  .  .  U  )  = 

j|  1  n' 


P(A  )P(Uj  A  )P(U2  A  ) 


P(U  A.) 
n  r 


The  terms  P(Aj)  and  P(Ui)  through  P(Uj^)  refer  respectively  to  the  prior  prob¬ 
abilities  of  the  target  class  and  of  the  areas  parameter  set  U .  .  .  U^.  Similar 
probabilities  for  the  presence  of  the  other  target  classes  within  a  given  area  will 
be  formed  from  the  set  of  parameters  associated  with  a  given  area.  Comparison 
of  the  values  for  the  several  target  probabilities  thus  obtained  will  indicate  to 
which  target  class  the  area  belongs.  When  the  probabilities  are  converted  to 
logs  to  facilitate  the  computations  and  the  prior  parameter  probabilities,  P(Ui) 
through  P(Un)i  are  omitted  (since  they  are  common  to  all  the  target  class 
probabilities),  a  target  recognition  vector,  Rj,  for  each  target  class,  Aj,  can 
be  formed  from  equations  50  and  51; 


i  =  n 

R.  =  log  P(A.)  -  >  A..  +  B..U.  +  C..U.^  .  (52) 

J  ®  '  j'  ij  ij  1  ij  1 

i  =  1 


A  similar  recognition  vector  will  be  formed  for  the  nontarget  background.  The 
largest  recognition  vector,  among  the  several  that  are  formed  from  a  given 
area,  will  identify  the  target  or  nontarget  class  to  which  the  area  most  probably 
belongs.  The  difference  between  the  two  largest  target  class  vectors  (since 
they  are  logs)  will  be  equivalent  to  the  ratio  of  the  corresponding  two  target 
class  probabilities  (equation  51)  and  therefore  will  be  proportional  to  the  odds 
that  the  given  area  belongs  to  a  given  target  class  rather  than  to  its  nearest 
competitor. 


Recognition  Vector  Profile 

The  photo-interpreter  for  this  program  is  required  to  locate  and  identify  one 
or  more  ♦^argets  in  a  given  picture,  so  that  in  effect  he  will  examine  each 
target-s.  .  subarea  for  the  probability  that  it  contains  one  of  the  targets  he 
knows  is  present.  Therefore,  a  model  that  is  to  predict  the  probability  that 
the  photo-interpreter  has  correctly  identified  the  targets  that  are  known  to  be 
present  must  evaluate  at  each  subarea  the  recognition  vectors  for  the  above 
known-to-be-present  target  classes.  The  sequence  of  recognition  vectors 
evaluated  at  each  subarea  for  a  given  target  class  form  ?  recognition  vector 
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profile  of  the  given  target  class.  The  recognition  vectors  of  a  given  target 
class  profile  are  compared  with  one  another  so  that  the  area  having  the  larges 
recognition  vector  most  probably  contains  the  target.  Even  though  the  para¬ 
meter  values  differ  from  one  area  to  the  next,  the  prior  parameter  probabili¬ 
ties  may  still  be  neglected,  since  these  probabilities,  P(Ui),  are  assumed  to 
be  constant  and  independent  of  the  value  of  Ui  obtained. 

The  difference  between  the  two  largest  recognition  vectors  of  the  recognition 
profile  is  proportional  to  the  odds  that  the  target  is  contained  in  a  given  sub- 
area  rather  than  in  its  nearest  competitor.  A  difference  can  also  be  found  fctr 
tiie  recognition  vector,  (Rt)i  that  corresponds  to  the  area  known  to  contain  the 
target,  and  for  the  largest  vector  among  all  the  other  subarea  vectors.  This 
latter  difference  represents  the  odds  that  the  target  is  indeed  at  the  target  area 
rather  than  at  its  nearest  competitor,  and  is  therefore  related  to  the  probability 
that  the  photo-interpreter  has  correctly  identified  the  given  target.  Other  use¬ 
ful  information  may  be  extracted  from  the  recognition  profile.  For  example, 
the  difference  between  the  target  area  vector  and  the  average  value  of  the  other 
vectors  is  a  measure  of  the  difficulty  the  photo-interpreter  has  in  locating  and 
identifying  the  given  target. 

Evaluation  of  the  Target  Detection  Probability 

The  actual  evaluation  of  the  target  detection  probability  (the  probability  that 
the  photo-interpreter  is  correct)  must  relate  the  varying  abilities  of  the  photo¬ 
interpreter  to  detect  the  targets  with  the  recognition  vector  results  obtained  by 
the  model. 


One  approach  is  to  assume  that  the  model-acquired,  target-class,  recognition- 
vector  differences  (Rtj  -  Rij)  differ  by  a  constant  (C)  from  the  equivalent  results 
obtained  by  an  "average"  photo-interpreter.  It  is  further  assumed  that  a 
Gaussian  distribution  function  describes  the  variations  in  the  target  detection 
abilities  of  the  various  photo-interpreters.  Assume  (1)  that  Qj;  represents  the 
difference  in  the  recognition  vectors  for  target  class  j  obtainea  by  a  photo¬ 
interpreter  at  the  target  subarea  and  at  a  given  no-target  subarea,  i  (a  correct 
decision  will  be  made  by  a  photo -interpreter  at  subarea  i  if  the  previously  de¬ 
termined  target  vector  is  greater  than  the  area  i  recognition  vector  or  ajj  >  0), 
and  (2)  that  Rj^j  -  Rjj  is  the  model-acquired  difference  between  the  recognition 
vectors  for  target  class  j  evaluated  (equation  52)  at  target  area  (Rtj)  and  sub- 
area  i  (Rij)-  The  probability  that  a  correct  photo-interpreter  decision  is  made 
with  respect  to  subarea  i  will  be 


This  probability  may  be  converted  into  the  simpler  error  function  form  if 


Then 
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dx  . 


P{a..  > 


This  in  turn  reduces  to 


and  the  values  of  the  Erf  function  may  be  found  from  tables  (for  example, 
reference  14,  p  116). 


Finally  the  target  detection  probability,  Pcj»  for  target  class  j  is  the  probability 
that  the  photo-interpreter  has  made  the  correct  no-target  decision  at  each  of 
the  N  no-target  subareas: 

i  =  N 

^cj  = 

i  =  1 


or 


(56) 


-log  P 


cj 


L  =  N 

log  P(a.^  >  0)  . 

L  =  1 


The  constants  C  and  a  associated  with  the  Erf  function  (equation  55)  reflect 
the  different  abilities  of  the  photo-interpreters  who  are  to  participate  in  the 
actual  psychophysical  experiments  and  are  to  be  obtained  from  the  results  of 
the  experiments. 


Parameter  Measurements 


To  ensure  that  the  parameter  set  effectively  represents  the  photographic  infor¬ 
mation  that  is  available  to  the  photo-interpreter,  as  many  parameters  as  pos¬ 
sible  should  be  measured.  The  parameters  that  can  be  extracted  from  a  given 
readout  signature  are  listed  under  c  above;  signatures  are  suggested  under  b  to 
describe  the  two-dimensional  photographic  characteristics  of  pattern,  texture. 
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and  intensity,  which  will  be  used  by  the  photo-interpreter  to  achieve  target 
recognition.  Circuitry  must  be  developed  to  obtain  the  suggested  parameters 
and  signatures. 

Metric  I  required  two  signatures,  each  of  which  was  obtained  in  the  x  and  y 
readout  directions:  the  photographic  transmissivity,  tj^  and  ty,  and  the  ab¬ 
solute  gradient  (the  absolute  value  of  the  first  derivative  oi  the  transmissivity). 

The  following  parameters  can  be  made  available  from  the  Metric  I  signature 
set  for  Metric  II  (reference  15,  PP  3-25): 


1 . 


The  average  transmissivities. 


<.t  > 


X 


and  <  t  >  , 

y 


2. 

3. 

4. 


2  2 

The  transmissivity  second  moments,  <t  >  and  <t  >  , 

X  y 


The  average  gradients,  < 


8t 


>  and< 


at 


>,  and 


The  gradient  second  moments,  < 


at 


'>and<||i|'>. 

ay 


The  first  and  second  moments,  since  the  sampling  interval  and  therefore  the 
number  of  sample  points,  N,  are  constant,  are  equivalent  to  the  signal  d-c 
content  and  energy. 


Two  signatures  involving  the  detection  of  target  edges  must  be  obtained  for 
Metric  II.  These,  obtained  in  both  the  X  and  Y  directions,  are  the  absolute 
gradients,  |G|,  that  exist  in  the  vicinity  of  the  target  edges,  and  the  widths, 

X,  between  the  successive  edges.  The  readout  video  will  be  restricted  so  that 
edges  are  obtained  only  at  the  very  high  or  very  low  transmissivity  regions  of 
the  photograph;  these  restrictions  are  applied  because  the  targets  to  be  con¬ 
sidered  will  have  either  a  very  high  or  very  low  transmissivity.  These  two 
signatures  provide  the  following  parameters; 


1.  Total  number  of  pattern  edges,  N-(x)  and  N  (y)  , 

6  o 


1  =  Ng(x) 


i  =  Ng(y) 


2.  D-c  content  of  gradient,  ^  ^  |G.(x)j  and 

i  =  1  i  =  1 


i  =  Ng(x)  i  =  N^(y)  2 

3.  Energy  of  gradient,  ^ 

i  =  1  i  =  1 


i  =  N. 


i  =  N 


X  y 

4.  D-c  content  of  pattern  widths,  (X.)  and  (Y.),  and 

i  =  1  i  =  1 
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5.  Energy  of  pattern  widths. 


1  =  N  I  =  N 

i  =  1  i  =  1 

These  parameters  for  Metrics  I  and  II  are  far  from  being  a  complete  repre¬ 
sentation  of  the  photographic  information  that  is  available  to  the  photo-inter¬ 
preter  because  they  contain  no  information  of  the  target  orientation  (only  data 
from  two  orientation  directions  are  obtained)  or  of  the  phase  or  interdependence 
of  adjacent  signature  measurements.  Other  signatures  have  been  suggested, 
such  as  edge  length  and  spacing  signatures  obtained  from  edges  that  have  a 
restricted  separation. 

Application  of  Recognition  Vector  Technique  to  Automatic  Target  Recogni- 
tion 


In  accordance  with  the  described  model,  the  probability  that  the  photo-inter¬ 
preter  correctly  identifies  the  target  within  a  given  picture  will  be  a  function 
(equations  54  through  56)  of  the  recognition  vectors  evaluated  at  the  various 
subareas  of  the  picture.  .An  experimental  program  involving  actual  perform¬ 
ance  of  the  photo-interpretation  task  would  test  the  validity  of  this  concept. 

Once  the  concept  of  the  recognition  vector  is  proved  to  be  valid,  the  application 
of  the  vector  can  very  readily  be  extended  to  actual  automatic  target-recognition 
decisions  that  will  alleviate  the  target  identification  tasks  of  the  photo-interpreter. 
The  actual  computation  of  the  recognition  vector  (equation  52)  from  a  set  of 
signature  parameters  can  be  accomplished  easily  since  only  linear  operations 
are  involved.  The  three  constants  (Ajj,  Bjj,  and  Cjj)  of  the  distribution  func¬ 
tion  for  each  target  class  and  for  the  nontarget  background  class  as  well  as 
for  each  parameter  can  be  stored  in  the  recognition  equipment  as  three  re¬ 
sistors.  Summation  of  the  parameter  signal  currents,  Uj  ,  applied  to  t!ie 
proper  resistors  will  provide  an  indication  of  the  target-class  recognitiori  vec¬ 
tor  for  a  given  subarea.  Electronic  comparison  of  the  various  target-class 
vectors  for  a  given  area,  to  determine  the  largest  vector,  will  permit  auto¬ 
matic  identification  of  the  area.  The  simplicity  of  the  parameter  measurement 
and  decisioning  process  to  form  a  recognition  decision  allows  the  entire  recog¬ 
nition  process  to  become  airborne. 

Summary 


A  model  has  been  described  that  utilizes  various  raster  readout  signatures  of 
the  individual  subareas  of  a  given  picture  matrix  to  evaluate  a  target  detection 
probability;  that  is,  the  probability  that  the  photo-interpreter  will  correctly 
identify  the  target  subarea. 

The  experimentally  acquired  signatures  measure  those  characteristics  as 
pattern  texture  and  intensity  that  a  photo-interpreter  ordinarily  uses  to  identify 
a  target.  Each  signature  of  a  given  subarea,  since  it  is  in  effect  a  sampled 
electrical  signal,  is  described  by  a  set  of  parameters  that  reflect  the  signal's 
electrical  characteristics  as  signal  length,  energy,  d-c  content,  and  spectral 
energy  distribution 

Statistical  distribution  functions  of  these  parameters  are  obtained  from  pre¬ 
viously  acquired  experimental  data  for  each  target  class  that  is  oT  interest  to 
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the  photo-interpreter.  Insertion  of  the  parameters  obtained  from  a  given  un¬ 
known  subarea  into  the  corresponding  parameter -distribution  functions  of  a 
given  target  class  results  in  the  formation  of  a  target -class  recognition  vector 
for  the  subarea.  This  vector  actually  is  the  log  of  the  probability  that  the  sub- 
area  contains  the  target  class.  Comparison  of  the  several  target-class  and 
nontarget  recognition  vectors  that  are  formed  at  the  given  subarea  from  its 
set  of  signature  parameters  permits  an  identification  of  the  area.  Because  of 
the  equivalence  between  this  identification  decision  and  that  rendered  by  the 
photo -interpreter  at  each  subarea,  the  recognition  vectors  obtained  at  all  the 
subareas  may  be  combined  (equations  54  through  56)  to  provide  the  desired 
target  detection  probability. 
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APPENDIX  U 


METRIC  m,  TARGET -IDENTIFICATION  TIME 
UTILIZING  STATIC  ANALYZER  MEASUREMENTS 

INTRODUCTION 

Target  recognition  in  Metric  III  was  treated  as  the  process  of  forming  a  suf¬ 
ficient  correlation  between  the  received  data  and  some  stored  or  a  priori  data 
possessed  by  the  observer.  A  description  of  the  recognition  process  must 
therefore  be  a  function  of  the  stored  data  as  well  as  the  received  data.  A  de¬ 
scription  of  data  stored  within  a  human  observer  is  virtually  impossible  to 
obtain.  The  next  best  thing  is  to  use  the  physical  entity  that  gives  rise  to  the 
stored  data  This  will  be  referred  to  as  the  reference  graphic.  This  essen¬ 
tially  is  a  pictorial  prediction  of  selected  portions  of  the  area  under  considera¬ 
tion.  At  some  point  in  the  reference  graphic,  a  target  is  located  and  the  ob¬ 
server  is  to  find  the  corresponding  point  in  the  complete  picture  --  the  problem 
scene.  Presumably,  if  the  reference  and  the  problem  scene  are  identical,  then 
the  observer  should  experience  the  minimum  difficulty  in  locating  the  target. 
However,  the  reference  may  differ  considerably  from  the  problem  scene  both 
in  completeness  and  accuracy.  Therefore,  the  difficulty  of  the  recognition  task 
appears  to  be  a  function  of  (1)  content  of  reference,  (2)  difference  between  ref¬ 
erence  and  scene,  and  (3)  quality  and  complexity  of  problem  scene.  Items  1  and 
2  are  not  independent,  since  the  quality  of  the  reference  could,  for  example,  be 
enhanced  at  the  expense  of  creating  a  large  difference  between  it  and  the  prob¬ 
lem  scene.  This  model  attempts  to  relate  these  ideas  to  the  time  an  observer 
requires  to  locate  a  target.  The  measurements  to  be  made  utilize  the  correla¬ 
tion  functions  of  the  reference  scene  and  the  problem  scene. 

DEFINITION  OF  MODEL  EQUATIONS 

Figure  19  represents  the  reference  wherein  the  origin  of  the  coordinate  system 
x',  y'  IS  placed  for  convenience  at  the  target  that  is  to  be  found  by  the  observer. 
Various  features  or  cues  are  indicated  in  the  reference,  and  the  target  may  or 
may  not  be  represented  by  a  feature.  The  problem  then  is  to  locate  the  cor¬ 
responding  point  in  the  scene  (B  in  figure  19);  that  is,  the  observer  is  to  find 
the  coordinates  x,  y  of  the  target  relative  to  an  arbitrary  coordinate  system 
in  the  problem  scene.  If  the  problem  scene  and  the  reference  did  not  differ  in 
content,  then  the  relationship  between  them  could  be  expressed  by 

s(x,  y)  =  r(x',  y') 

=  r(x  -  u,  y  -  v)  ,  (57) 

where  x  =  u,  y  =  v  are  the  coordinates  of  the  target  in  the  problem  scene. 
Equation  57  states  that  the  scene  and  the  reference  look  alike  except  for  a 
shift  u,  V  between  their  respective  coordinate  systems.  In  general,  the  scene 
and  the  reference  will  not  have  the  same  content,  so  that  in  addition  to  a  co¬ 
ordinate  shift  there  will  be  a  difference,  n(x,  y),  that  must  be  added  to  the 
referent..e  before  the  two  would  be  equal: 

s(x,  y)  =  r(x  -  u,  y  -  v)  +  n(x,  y)  .  (58) 

The  function  n(x,  y)  contains  all  the  differences  in  content  between  the  problem 
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scene  and  the  references  and  is  actually  defined  by  equation  58;  thus, 

n(x,  y)  =  s(x,  y)  -  r(x  -  u,  y  -  v)  .  (59) 

The  problem  of  finding  the  target  may  now  be  stated  in  terms  of  solving  equa¬ 
tion  56  for  u  and  v.  Thus,  the  observer,  having  been  briefed  with  r(x',  y'), 
now  inspects  s(x,  y),  seeking  satisfactory  values  for  u  and  v.  Since  the  term 
n(x,  y)  is  not  known,  the  only  recourse  is  to  give  it  a  probabilistic  interpreta¬ 
tion.  This  is  equivalent  to  the  statement  that,  although  n(x,  y)  is  unknown, 
certain  constraints  upon  it  are  assumed.  Solutions  for  u  and  v  therefore  will 
be  probabilistic.  The  target  coordinates  x  =  u,  y  =  v  can  never  be  known  with 
certainty.  It  is  only  possible  to  maximize  the  probability  that  the  correct  values 
have  been  selected. 

Dependent  upon  the  assumptions,  a  probability  distribution  for  the  coordinates 
of  the  target  is  obtained.  Associated  with  this  final  distribution  is  a  gain  of 
information  due  to  the  process.  It  is  proposed  that  the  model  relate  the  recog¬ 
nition  time,  t,  to  the  gain  of  information,  Ig;  thus, 

t  oc  ^  .  (^>0) 

g 

In  addition  to  the  information  gain,  aterm,  C(s),  expressing  the  complexity  of 
the  scene,  is  included.  This  is  a  measurement  of  the  detail  in  the  scene.  The 
expression  for  t  therefore  becomes 

t  =  PJ-+  C(s)  .  (61) 

g 

Two  cases  are  considered,  deper.ding  upon  the  form  of  the  reference  material. 
These  are  continuous -level  and  three-level  (black,  gray,  and  white).  Equa¬ 
tion  61  is  developed  as  follows  for  each  of  these  cases. 


□ 

0  0. 

-  -  ' 

■ 

target  ^ J 

w 

E 

4 

Figure  19  -  Representation  of  Metric  III  Target  Reference 
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CONTINUOUS-1.EVEL  REFERENCE 


The  variance  of  u  is  given  by 


and  for  v  by 


where 


- ; - - rr/i 
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The  information 


=  ^(0,  0)  =  the  amplitude  of  the  peak  of  the  autocorre- 
lated  reference, 

=  the  second  derivative  of  the  autocor related  reference 
with  respect  to  the  u  direction  evaluated  at  the  peak, 

=  the  second  derivative  of  the  autocorrelated  reference 
with  respect  to  the  v  direction  evaluated  at  the  peak, 

=  the  area  of  the  reference, 

=  the  reference  to  noise  power  in  the  scene,  and 

=  a  constant. 

gain  represented  by  equations  6Z  and  63  is  given  by 


I 

g 


A  R.^  "(0,  0)^  "(0,  0) 
,  r  1  u  '  '  ' 

log - 2 - 

[k  Mo,  0)] 


(64) 


The  complexity  of  the  problem  scc>ne,  C(s),  is  measured  by  the  logarithm  of 
its  normalized  derivative  energy. 


C(s)  =  o)[>^u"<0.  0)  +  0)|  , 


(65) 


where  /  denotes  the  problem  scene  autocorrelation  fui.ction.  The  model  for 
this  case,  equation  61,  therefore  becomes 


t  =  k  log 


k[/(0,  0)] 


(0,  0)/i^"  (0,  0) 


log 


>)(0,  0) 


K"  '«• 
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.  (0, 
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(66) 
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THREE-1.EVEL  REFERENCE 


The  standard  deviation  of  u  is  given  by 


a  u 
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where 


^^■'(0,  0) 


w 

X 


w 

y 
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the  magnitude  of  the  slope  of  the  autocor related  reference 
taken  in  the  u  direction  evaluated  at  u  =  0  +  , 

the  minimum  width  of  cues  in  the  x  direction, 
the  minimum  width  of  cues  in  the  y  direction,  and 
a  constant. 


For  the  v  direction,  the  standard  deviation  is  given  by 


a  V 
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w  w 
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(68) 


where  (0,  0)  is  the  magnitude  of  the  slope  of  the  autocorrelated  reference 
taken  in  the  v  direction  evaluated  at  v  =  0+. 


The  information  gain  represented  by 
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log  ' 
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equations  67  and  68  is  given  by 
(0.  0)^^'  (0,  0) 
k[^(0,  0)] 


(69) 


Again,  by  use  of  equation  66  for  the  complexity  of  the  scene,  the  average  time, 
from  equation  6Z,  becomes 
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t  =  k^  log  ^Kf  ^  1  ^  f  k,  log 
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The  undetermined  constants  in  equations  66  and  7  0  would  be  obtained  by  a 
least-squares  fitting  process  similar  to  that  discussed  for  Metric  I. 
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APPENDIX  m 


ANALOG  COMPUTER  CIRCUITRY 


INTRODUCTION 

The  four  transmissivity  characteristics  of  the  pictures  selected  for  Metric  1 
were  measured  by  a  flying  spot  scanner -analog  computer  setup.  The  com¬ 
puter  programmed  the  scan  of  the  flying  spot  and  read  out  the  measurements 
of  each  characteristic  on  an  individual  channel.  Each  picture  was  divided  into 
256  equal  areas,  16  vertical  and  16  horizontal.  Each  area,  or  cell,  was 
scanned  10  times  horizontally  and  10  times  vertically  (fast  scan).  The  cells 
were  scanned  column  by  column  and  row  by  row  (slow  scan). 

SCAN  CONTROL  CIRCUITRY 

The  fast  scan  is  the  basic  clock  pulse  train  that  times  and  synchronizes  the 
scan  control.  It  is  generated  by  means  of  a  group  of  computer  amplifiers  (18, 

19,  and  20)  in  figure  20.  Because  of  the  limiter  branches  in  its  feedback,  the 
output  of  amplifier  20  will  be  either  +20  or  -20  v.  This  output  is  compar-'d  to 
the  ramp  output  of  amplifier  19  at  the  input  of  amplifier  18.  Each  time  the 
ramp  reaches  the  20-v  level,  the  output  of  amplifier  18  changes  polarity.  This 
flips  the  output  of  amplifier  20  to  the  opposite  state.  The  gains  of  the  two  input 
branches  into  amplifier  19  are  adjusted  to  give  a  slowlv  rising  ramp  of  one 
polarity  and  a  fast-climbing  ramp  of  the  other  polarity,  the  last  corresponding 
to  the  fly-back  time  between  scan  lines.  The  upper  trace  of  figure  21  shows 
the  wave  shape  involved. 

The  short  pulse,  as  reflected  at  the  output  of  amplifier  20,  is  used  to  generate 
the  slow-scan  by  means  of  integrator  24  and  comparator  22.  When  the  stair¬ 
case  output  of  24  reaches  the  voltage  set  by  potentiometer  13,  relay  R1  is  ener¬ 
gized  and  momentarily  discharges  the  integrator  back  to  zero.  This  happens 
every  11  steps,  as  shown  by  the  middle  tracing  of  figure  21.  Two  parallel  blank¬ 
ing  circuits  energized  by  the  fast  and  slow  scans  generate  a  short  and  a  long 
blanking  pulse  which,  after  being  shaped  and  properly  biased  in  amplifiers  10 
and  11,  are  applied  either  to  the  cathode  of  the  scanner  by  amplifier  11  or  to 
electronic  switches  at  the  computer  operating  on  the  video  signal  (amplifiers  4 
and  4a).  The  short  pulse  exists  at  the  beginning  and  at  the  end  of  each  fast  scan. 
The  long  pulse  covers  the  transient  time  between  one  slow  scan  and  the  next, 
obliterating  one  of  the  11  steps.  Thus,  each  elemental  square  of  the  picture  is 
covered  by  10  scan  lines  in  each  direction.  Blanking  is  applied  at  the  scanner  if 
the  end  product  is  to  be  a  picture  or  at  the  electronic  switch  if  the  end  product 
is  to  be  a  voltage.  The  blanking  pulses  are  shown  at  the  lower  graph  of  figure  21. 

The  blanking  circuits  (figure  22)  combine  a  level  detector  and  a  timed-pulse 
generator.  They  were  biiilt  as  separate  units  and  added  to  the  computer.  The 
front  part  of  the  circuit  is  a  Schmitt  trigger,  the  level  of  which  can  be  adjusted 
by  means  of  a  potentiometer.  This  establishes  the  timing  of  the  leading  edge 
of  the  blanking  pulse.  The  output  of  the  level  detector  (trigger)  is  applied  to 
the  adjustable -delay  multivibrator  shown  at  the  right-hand  side  of  the  diagram. 
There,  the  duration  of  the  generated  pulse  is  controlled  by  capacitor  Cd.  The 
durations  of  intraline  and  interelement  blankings  are  established  by  properly 
selecting  the  values  of  these  capacitors. 

The  same  pulse  from  relay  R1  that  discharges  slow-scan  integrator  24  also 


triggers  a  bistable  multivibrator  circuit,  called  the  alternator,  whose  alternat¬ 
ing  output  pulses  energize  relays  RZ  and  R3.  These  relays  operate  indicator 
lights,  R  and  Y,  and  a  counter  C,  and  also  switch  the  fast  and  slow  scans  alter¬ 
nately  between  the  X  and  Y  deflection  coils  of  the  scanner  via  amplifiers  11  and 
12.  Through  the  same  amplifiers,  constant  bias  voltages  are  added  (potenti¬ 
ometers  4  and  5)  to  position  the  spot  properly  at  one  of  the  corners  of  the  pic¬ 
ture;  also,  either  the  horizontal  or  the  vertical  steps  are  introduced  here. 
These  steps  originate  at  the  alternator,  which  has  two  plate  voltages.  One 
voltage  is  differentiated  to  produce  two  pulses  of  opposite  polarities  per  switch¬ 
ing.  One  pulse  is  isolated  through  a  biased  diode  and  fed  into  the  X-step  gen¬ 
erator,  amplifiers  7  and  15.  These  amplifiers  are  arranged  in  the  same 
manner  as  amplifiers  22  and  24;  hence  a  staircase  output  is  produced  by  am¬ 
plifier  1  1 . 

The  comparator  amplifier,  15,  is  set  through  potentiometer  14  to  operate 
when  the  alternator  has  completed  16  periods  (when  the  16  square  elements 
of  one  row  on  the  picture  have  been  scanned).  At  the  end  of  this  interval, 
relay  R4  discharges  integrator  7  and  at  the  same  time  increases  the  output 
of  integrator  21  by  one  step,  thus  moving  the  scan  one  cc’uma  over.  As  shown 
in  figure  20,  this  output  is  added  to  the  input  of  the  other  deflection  amplifier, 
12.  Comparator  amplifier  24B,  in  series  with  relay  R4,  is  adjusted  to  en¬ 
ergize  relay  R5  at  the  end  of  16  steps  out  of  integrator  21.  This  throws  the 
computer  to  the  "hold"  position,  freezing  all  voltages  and  terminating  the 
process. 

Figure  23  shows  ii  succession  the  slow  scan  produced  by  amplifier  24;  the 
pulse  from  relay  R1  into  the  alternator;  one  of  the  alternator  outputs;  the 
pulse  that  is  part  of  the  derivative  of  this  output  and  feeds  integrator  7;  the 
output  of  this  integrator  (horizontal  steps);  and  the  output  of  integrator  21 
(vertical  steps). 

Figure  24  shows  the  alternator  circuitry.  The  upper  part  is  a  monostable 
multivibrator  actuated  by  the  pulse  from  relay  Rl.  Its  purpose  is  to  generate 
a  reliable  pulse  of  constant  length  undisturbed  by  the  noise  peaks  of  the  input 
pulse  arising  from  contact  noise.  The  output  of  this  part  triggers  the  bistable 
multivibrator  shown  in  the  lower  part  of  the  figure,  which  turns  on  and  off  rl- 
ternately  a  red  and  a  yellow  lamp  indicating  to  the  equipment  monitor  the  di¬ 
rection  of  the  fast  scan  (X  or  Y);  it  also  steps  counter  C,  which  indicates  the 
position  of  the  elemental  square  last  scanned. 

The  outcome  of  the  scan  control  circuitry  described  above  is  the  coverage  of 
a  square  picture  area  with  a  length,  L,  on  each  side  of 

L  =  2.  25  in.  (71) 

Z 

by  a  matrix  of  N  square  elements,  where 

N  =  16,  (72) 

each  having  side  length  of 

/  =  0.  140  in.  (73) 
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Figure  20  -  Analog  Computer  Scan- Control  Circuit 


Figure  22  -  Blanking  Circuit 
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Each  element  in  turn  is  scanned  by  a  number  of 


n  =  10  (74) 

equidistant  parallel  lines  in  the  X  direction  and  by  as  m.any  in  the  Y  direction. 
The  duration  of  each  scan  line  is 


tj  =  0.  0/2  sec  .  (75) 

and  therefore  the  writing  speed  at  the  picture  plane  is  very  nearly 

i 


=  2  in.  per  second  .  (76) 

A  blanking  interval  of  0.  024  sec  is  provided  between  adjacent  lines  and  0.  120 
sec  between  elements  and  between  rows. 

An  additional  feature  of  the  scan  control  circuitry  is  the  brightness  control 
loop  shown  at  the  lower  right-hand  part  of  figure  20.  In  this  loop,  the  output 
of  a  monitor  photocell  facing  the  cathode  ray  tube  is  continuously  compared  to 
a  reference  voltage  (potentiometer  P30)  at  the  input  of  amplifier  30B.  The 
output  of  this  amplifier,  connected  to  the  control  grid  of  the  crt,  keeps  the  light 
flux  emanating  from  the  flying  spot  at  a  constant  level. 


77 


SIGNAL  MEASURING  CIRCUITRY 


The  electrical  output  of  a  photomultiplier  tube  (the  video  signal,  v)  is  channeled 
to  the  computer  via  amplifier  lb,  which  is  contained  in  the  scanner,  and  a  gain- 
control  potentiometer,  P,  inserted  between  the  amplifier  output  and  ground. 

The  signal  is  pres  mted  to  the  computation  circuit  (figure  25)  by  two  buffer 
amplifiers,  14  and  17A,  which  have  t)\e  transfer  function 


-k 


1 


(77) 


with 


=  5  , 

=  3.  35  X  10'^  sec,  and  (78) 

4  =  0.5. 

The  amplifier  therefore  behaves  very  much  like  a  simple  low-pass  filter  with 
a  3-db  cut-off  point  at 


f ^  =  47  50  cps  , 


(79) 


corresponding  to  spatial  cut-off  frequency  at  the  picture  plane  of 

ru 

=  3165  cps  per  inch  , 


(80) 


where  u  is  given  by  equation  76  and  r  is  the  magnification  factor  of  the  optical 
path  between  picture  plane  and  crt  face;  that  is, 


r  =  ^-  =  0.  75 
o 


(81) 


The  above  frequencies  provide  a  margin  for  additional  filtering,  rendering  the 
circuit  less  sensitive  to  noise  and  stray  pickups.  This  is  so  because,  assuming 
a  practical  spot  size  of  zq  =  2.  5  mils  at  the  crt  face,  or 


z  =  rz  =  1 . 87  mils 
o 


(82) 


at  the  picture  plane,  the  cut-off  point  because  of  spot  size  will  be 


~  -  530  cps  per  inch 


(83) 
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Figure  25  -  Signal  Proceasing  and  Measuring  Circuit 


(84) 


=  Uf^  .  =  1,  060  CDS  . 


Thus,  the  crt  spot  contribution  may  be  described  by  the  approximate  transfer 
function  of 


(85) 


with 


^2  277f- 


=  15X  10  ^sec 


(86) 


The  additional  filtering  was  implemented  by  means  of  a  33-^fif  capacitor  across 
the  2 -meg  feedback  of  amplifier  14  or  17A,  corresponding  to  a  time  constant  of 
^3  =  6.  6  X  10"^  sec  and  a  cut-off  frequency  of  2410  cps,  or  about  1200  c/cles 
per  inch.  Hence,  equation  77  has  to  be  changed  to 


1 


T 


2  2 


(87) 


(Tj  s  +  2C  s  +  l)(r3S  +  1) 


with 


kj  =  5  , 


-5 


(88) 


Tj  =  3.  35  X  10  sec  , 

^  =  0.  5,  and 

_5 

=  6. 6  X  10  sec  . 

The  measured  response  characteristics  are  shown  in  figure  26. 

serves  to  derive  the  four  parameters  of  interest  for  the  present  effort, 
namely  the  mean  and  variance  of  transmissivity  and  of  detail,  all  of  which 
have  already  been  discussed. 

In  terms  of  voltages  over  each  elemental  square,  these  quantities  will  be 

h  . 

dt  ,  (89) 


tl 


m  -tj 

=  ^  Vi*  +  V,(x.,  y)| 

i  =  \  Jo 

m  .t 

E  /  1^1 

‘  i  =  1  •'o 


^  (x,  y^)  +  (x.,  y)  dt  - 


(90) 
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PHASE  SHIFT  DEGREES 


m 


^  1  r 

>‘di^  nh-u  ^  I  ''i 

i  =  1  •'o 


m  -t 


^di  j  2mt ,  u 


i  =  1  •'o 


v,(x,  yM  + 


h'. 


and  also 


a  =  =  0.  0196 


q'  =  2mtj  =  1.40  sec  , 


(91) 


(92) 


(93) 


where  a  is  the  area  of  the  elemental  square,  and  O' '  its  equivalent  in  time  units. 
Accordingly,  four  parallel  branches  are  needed  to  compute  s^.,  p  j--  3.nd  s^. 

concurrently. 


To  measure  p^.,  integrator  9  of  figure  25  is  used  after  provisions  are  made  to 


FREQUENCY  (CYCLES  PER  SECOND) 

Figure  26  -  Amplifier  Characteristics 
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blank  electronical .y  between  the  ending  of  one  line  and  the  beginning  of  the 
next,  as  well  as  between  square  elements.  I'o  this  end,  the  positive  blanking 
pulse  supplied  by  amplifier  4,  figure  20,  is  used  to  annihilate  the  always 
negative  Vj  at  the  input  of  the  integrator  amplifier  9  by  driving  to  ground  the 
diode  shown  in  series  with  this  amplifier  output.  A  0.2-v  bias  voltage  would 
have  been  necessary  at  the  same  output  to  compensate  for  the  diode  dead-space 
characteristic.  However,  a  similar  diode  across  the  feedback  path  of  the  am¬ 
plifier  causes  its  output  to  assume  either  a  +0.  2-v  or  -125-v  value. 

The  same  technique  is  employed  (via  amplifier  4a)  in  measuring  s^.^,  for  which 
Vj  is  squared  in  an  electronic  multiplier  prior  to  being  integrated  in  amplifier 
2.  At  the  same  time,  the  output  of  integrator  9  is  being  squared  and  subtracted 
from  that  of  integrator  2,  thus  yielding  «ti  dire  ctly.  The  frequency  and  phase 
characteristics  of  the  particular  multiplier  used  as  actually  measured  are 
shown  in  figure  27.  The  amplitude  and  phase  response  of  the  four -quadrant 
unit  to  the  product  of  a  100-v  d-c  and  a  sinusoid  voltage  of  25  v  peak-to-peak 
is  given  by  the  broken  lines.  The  solid  line  shows  the  response  of  one  quad¬ 
rant  to  a  half-wave  rectified  sinusoid  of  50-v  maximum  amplitude.  The  latter 
arrangement  is  the  manner  in  which  squaring  is  obtained  in  the  present  setup. 
The  d-c  characteristic  of  a  quadrant  is 


2 

with  <1  maximum  error  of  2  percent  of  the  measured  value,  y. 


Figure  27  -  Multiplier  Frequency  and  Phase  Characteristics 
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The  third  picture  parameter  ,  requires  the  differentiation  of  V ^  before 
averaging  over  each  elemental  square.  This  operation  is  performed  in  am¬ 
plifier  20A  which,  besides  the  unavoidable  input  time  delay  needed  in  all  elec¬ 
tronic  differentiators,  was  supplied  with  an  additional  delay  term  at  its  feed¬ 
back  to  attenuate  frequencies  higher  than  2400  cps.  The  resulting  transfer 
function  is 


where 


(T 


(94) 


k^  =  0.  0005  , 

=  0.  00008  sec,  and  (95) 

=  0. 000017  sec  . 


The  measured  gain  and  phase  characteristics  of  amplifiers  17A  and  20A  com¬ 
bined  are  given  in  figure  28.  On  the  other  hand,  as  shown  in  figure  29,  the 
3-db  cutoff  point  of  amplifier  lb,  which  follows  the  photomultiplier  in  the 
scanner,  occurs  at  a  very  high  frequency  for  this  amplifier  to  have  any  ef¬ 
fect  on  the  computations. 

By  contrast  to  and  s^j,  the  derivatives  aV^/dx  and  dVi/dy  may  be  either 
positive  or  negative.  Hence,  after  differentiation,  the  derivative  is  full -wave 
rectified  in  amplifiers  18A  and  19A  before  being  fed  into  integrator  1,  where 
electric  blanking  is  also  applied. 

The  second  moment  s  ..  of  detail  is  obtained  in  a  manner  similar  to  s^.  All 
four  integrators,  9,  1,  and  13,  are  momentarily  discharged  to  zero  ai  the 

end  of  the  complete  scan  of  each  elemental  square.  This  is  done  by  means  of 
relays  R6  and  R7  actuated  by  the  same  pulse  that  steps  up  the  vertical  -  step 
amplifier  7  through  amplifier  2a.  The  insertion  of  relays  R14  and  R15  allows 
for  some  time  delay,  which  is  desirable.  In  addition,  the  same  relays  keep 
the  integrator  outputs  to  zero  by  means  of  an  additional  input  to  amplifier  2a, 
effected  by  the  computer's  utility  switch  when  in  "hold"  position. 

As  shown  in  figure  25,  signals  A  and  B,  which  represent  V[  and  Vj  ,  are 
fed  into  amplifier  3  for  subtraction.  The  output  is 
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PHASE  SHIFT  (OECPEES) 


Figure  28  -  Differentiator  Characteristics 


Figure  29  -  Photomultiplier  Amplifier  Characteristics 
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where  T  is  the  effective  integration  time  per  elemental  square  when  blanking 
is  taken  into  account.  The  setting  of  the  gain  0  =  0.  1/R,  where  R  is  the  input 
resistor  associated  with  input  B,  is  done  on  the  iiasis  of  the  condition 

6^=0  when  =  K  =  constant  .  (9?) 


from  which 


1.  920 
T  • 


(98) 


Since  from  equations  74  and  84 

T  =  2mtj  =  1.44  sec  ,  (99) 

it  follows  that  R  should  have  the  value 

R  =  1. 333  .  (100) 

Of  course,  in  view  of  the  tolerances  in  all  the  condensers  and  resistors,  the 
exact  setting  of  R  should  be  done  experimentally  on  the  basis  of  equation  97. 

Similar  reasoning  is  used  in  the  case  of  amplifier  5a,  where  the  subtraction  of 

v/  ana 

is  done.  Figure  30  gives  a  sample  of  the  recorded  outputs  of  the  four  branches 
of  the  measuring  circuit. 

The  computations  shown  below  establish  the  values  of  the  mean  and  variance 
of  transmissivity  and  detail  over  a  cell,  and  of  the  mean  and  second  moment 
of  transmissivity  and  detail  over  the  whole  picture.  With  reference  to  fig¬ 
ure  25,  it  will  be  convenient  to  establish  a  factor  that  represents  the  input 
voltage  level  to  the  computing  circuitry  and  that  is  carried  through  all  compu¬ 
tations.  Thus  the  voltage  E'  out  of  the  control  potentiometer  P  is 

E'  =  - - - - ,  (101) 

1-  5  (1  -  i3) 

where  E  is  the  photomultiplier  output  and  the  potentiometer  setting.  On  the 
other  hand,  the  output  E  is 

E  =  0t  ,  (102) 

where  t  is  the  transmissivity  of  the  plate  at  the  location  of  the  flying  spot  and 
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0  =  volts  per  100  percent  transmissivity 


(103) 


is  the  photoelectric  gain  for  a  fixed  setting  of  the  crt  brightness  and  of  the 
brightness  control  feedback  loop. 

Hence, 


L'  =  - t  .  (104) 

^  -  1.  5(1  -  13) 

Throughout  the  measurement  period,  frequent  checks  of  0  were  made  by  in¬ 
serting  a  neutral  density  filter  with  a  transmissivity  of  39.  81  percent  at  the 
face  of  the  plate  and  measuring  E'  =  (E')^. 


amplifier  *a  <SRE  FIGURE  251 


Figure  30  -  Sample  Recording  of  Picture  Parameters 
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riie  r  fl'o  r  L“ 


1.5(1  -,i)  ,  3,  „ 

- 0 - — iw 

m 


This  quantity  is 


K 


h  t  1.  5(1  -  (1) 
0 


and  is  measured  in  terms  of  E'  as 


K  = 


0.  398 
(E') 

m 


With  this  in  mind,  the  following  relationships  are  found  (the  amplifier 
are  used  as  subscripts  to  their  output  voltages). 

1.  Average  Transmissivity  over  a  Cell: 

10 

‘'’6a  1 .  5  ^9  ’ 


10 

rrrr? 


4-  I  V 

:t  4  I  '^i 

•'o 


dt  , 


E  dt  . 


Since  ^ 


=  ‘“T 


0 


y  1. 5(1  -  /j) 


10 

'^ti  =  IT  ^i  =  ^'^i  ’ 


or 


=  0.  10  K(aWj)  transmissivity  . 


(105) 


(106) 


(107) 

numbers 


(108) 


Z.  Standard  Deviation  of  Transmissivity  over  a  Cell: 
From  equations  96  and  100, 


1 


'3a  (100)(11.  Ill) 


‘T  I 


-  —  /  V  ' 

T  J  M 

^  r\ 


dt 


(lOO)(l  1  111) 


-  m 


2V 


I 

>]• 


10.  000  v^* 


e-,  = 


3  0.  405  10,  000  Vj  ’ 


or 


=  =  '="'2  • 


=  90  bw^  ; 


since 


V,  =  5E'  = 


1 

50 


+  l.b(l  - 


E  . 


1 


+  1.  5(1  -  /3) 


t  = 


^1 


K  \  ^  Z 

Sy  =  3.  6K  (bw^)  transmissivity 


3.  Average  Detail  over  a  Cell; 


(109) 


10 

'23  1.5  ^1 


-  10 


-'O 


23 


=  2^ 


V. 


89 


0.  0005 


V.  -  |e 


20a| 

=  0. 0025 

=  0. 0025 

0. 0025 
K 


1 


1.  5(1  -  ,1) 

e 

4+  1-5(1  -13) 


dt  17A 


~  E 
dt 


=  0.  0025 


-Ie' 

dt  ^ 


dt 


But 


d  , 

d 

dt  ^ 

dx 

t 

^  dt 

—  u 

where  u  is  the  speed  of  the  flying  spot,  and  by  equation  76 

=  2d  . 


±t 

=  2 

-It 

dt 

dx 

,,  _  0,  005  ^ 

^2  "  -TT”^  ’ 

0.  00  5 


K  ^di 

0.  010 

II 

-o 

K  l^di 

t 

II 

100  K(.w 

3 


change  ner  inch  . 

4.  Standard  Deviation  of  Detail  over  a  Cell; 

By  analogy  with  e^^  in  2  above, 


5a  10,  000  ■  V. 


5  9  .  _  1 

5  0.  045  10,  000  10  ^^4’ 

iL  L 


0.  005  ^ 

d  , 


(110) 


9  0 


s 


0.  005^ 


=  4  X  lO^K^(fw)  . 

In  the  above  expressions,  a,  b,  c,  and  f  are  the  settings  of  the  recorder  chan¬ 
nels  in  volts  per  millimeter,  and  w^,  W2,  W3,  and  are  the  corresponding 
readings  in  millimeters. 
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APPENDIX  IV 


INFORMATION  CONTENT  OF  METRIC  U  TARGET -IDENTIFICATION 
PARAMETERS  AND  EVALUATION  OF  THEIR 
PROBABILITY  DENSITY  FUNCTION 


INTRODUCTION 

The  development  of  Metric  II  (see  appendix  I)  to  obtain  the  probability  of  cor¬ 
rect  tar  get -identification  required  quantative  expressions  of  the  information 
content  and  of  the  statistical  distribution  of  the  recognition  parameters.  These 
expressions  were  developed  as  explained  below. 

INFORMATION  CONTENT  OF  SIGNATURE  PARAMETERS  IN  RECOGNITION 
DECISION 

Derivation  of  Information -Content  Expression 


A  quantitative  expression  of  the  information  contributed  by  a  parameter  to  a 
recognition  decision  is  obtained  heur istically  and  corroborated  experimentally 
in  reference  16.  The  following  analysis,  based  upon  information  theory, 
yields  the  same  expression.  However,  the  expression  is  modified  to  obtain 
a  clearer  insight  into  the  role  that  a  parameter  plays  in  the  making  of  a  recog¬ 
nition  decision.  The  amount  of  information,  I(Ujj^),  delivered  by  a  particular 
value,  Ujj^,  of  a  signature  parameter,  U-,  to  the  recognition  of  a  target  c  lass, 
Aj,  is  given  by  (reference  17,  p  4) 

P(A  U  ) 

=  '08— 

J 


or  by  an  equivalent  expression  obtained  from  Bayes'  Rule  (  reference  1^, 
p  318). 


1 K 


(112) 


with 

P(A.)  =  the  a  priori  probability  that  a  particular  pattern  is  of 

class  A  . , 

J 

P(A^  ~  conditional  probability  that  the  unknown  pattern  is 

of  class  Aj  after  a  particular  parameter  measurement, 
U^j^,  of  the  pattern  has  been  made  available, 

P(U.j^)  =  the  probability  that  the  given  parameter  measurement 
is  obtained,  and 


P(Ui^  A^) 


the  probability  that  the  given  parameter  measurement 
is  obtained  from  a  particular  target.  A.. 
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If  P(Aj,  U^j^)  represents  the  distribution  of  the  ta  rget -a  ssociated  parameter 
values  over  the  entire  sample  of  N;  target  classes  and  Nj^  discrete  parameter 
values,  the  average  information,  itUj),  delivered  by  all  of  these  parameter 
values  becomes,  from  equation  112, 


I(U^) 


P(A^.  log 


(113) 


This  expression,  similar  to  that  obtained  in  reference  16,  is  modified  as  fol¬ 
lows.  The  log  ratio  is  converted  to  a  log  difference  and  the  following  substitu¬ 
tion  is  made. 


The  information  content  of  the  parameter  then  becomes 


P(A^)P(Uy^|A^)  log  P(U,^|Aj)  - 

k  = 

^  P(U.k)  ‘og  P(Vk) 

k  =  1 

In  terms  of  the  familiar  information  theory  entropy  function, 

H(x)  =  -  ^P(x)  log  P(x)  .  (114) 

X 

The  final  expression  for  the  information  content  of  the  parameter  (U.)  is 

j  =  N 

I(Ui)  =  -  ^  [p(A^)H(U^|A^)j  +  H(U^)  .  (115) 

J  =  1 

Thus,  the  average  information -contributing  capability  of  a  parameter  is  ex¬ 
pressed  by  the  difference  between  the  entropy  or  uncertainty  of  the  complete 
ac  ross -target  parameter  distribution  and  the  average  entropy  of  the  individual 
inuratarget  distributions  of  the  parameter.  Such  a  definition  agrees  with  in¬ 
tuitive  notions  for  desirable  target  recognition;  that  the  individual  intratarget 
distributions  have  as  low  an  entropy  as  possible  ( zero -entropy,  single -valued 
intratarget  distribution  functions  are  ideal)  and  that  the  over-all  distribution 
functions  of  the  parameter  over  all  the  intratarget  functions  be  as  broad  as 
possible  (have  high  entropy)  to  permit  discrimination  between  targets. 

Information  Contribution  of  a  Computed  Parameter 


J  =  N,  Nk 

E  E 

j  =  1  k  =  1 


In  some  cases  two  individual  parameters  may  be  combined  to  form  a  third 
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parameter.  For  example,  the  average  power  (Q)  of  a  sampled  signal  is  the 
quotient  of  the  energy  (E)  and  the  number  (S)  of  sampled  points.  If  the  para¬ 
meters  E  and  S  are  used,  what  additional  information  is  available  from  the 
third  computer  parameter,  Q?  In  a  given  pattern,  the  two  measured  para¬ 
meter  values  (Ej  and  S;)  determine  the  value  for  a  third  parameter  (Qj^)  The 
resulting  distribution  function  expressing  this  dependency  becomes  a  delta 
function: 


P(Q, 


E, 

1,  when  s-  =  “k 


0,  when  ^  ^  . 

j 
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The  resulting  entropy  therefore  is 


H(Q  E  S  ) 
1  J 


E^Sj)  log  P(Q^ 


E.S.)  =  0  . 
1  J 


(117) 


The  joint  entropy  of  all  three  parameters  (from  -eference  18,  p  5)  is 


H(QES)  =  H(ES)  + 


i  J 


E.S.) 
1  J 


which,  from  equation  117,  permits 


(118) 


H(QES)  =  H(ES)  .  (119) 

In  a  similar  manner  for  every  target  class, 

H(QES|Ay  =  H(ES|A^)  .  (120) 

Therefore,  from  equation  115,  the  information  content  of  all  three  parameters 
is  equal  to  that  of  two  alone,  or 

I(QES)  =  I(ES)  ,  (121) 


and  the  third  formed  parameter,  Q,  provides  no  additional  information. 


EVALUATION  OF  PARAMETER  PROBABILITY-DENSITY  FUNCTION 


Consideration  of  Single -Cluster  Data 

Experimental  values  of  the  probability  density  distribution,  f(u),  of  a  given 
parameter,  u,  for  <  u  <  are  available  for  a  given  target  class. 

The  distribution  of  these  values  is  normalized  as  shown  in  equation  122.  It 
is  required  that  this  data  distribution,  f(u),  be  represented  by  a  probability 
density  function  as  shown  in  equation  123. 
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with 


u  =  u 

max 

Y,  '<“>  =  ‘ 

11  =  u 

m  in 


P(u) 


+  bu  f  cu  ) 


u  =  u 

max 

E 

u  =  u 

min 


P(u) 


1 


(122) 


(123) 


(124) 


In  addition,  the  interval  between  the  various  sampled  values  of  u  must  be  con¬ 
stant  not  only  within  a  given  target  class  function  but  also  among  all  the  target 
class  functions,  so  that  tlie  various  values  of  P(u)  for  the  different  values  of  u 
and  for  the  different  target  classes  are  compatible  with  one  another. 

Consider  the  following  relationship,  D(u),  between  the  experimental  set  of 
data,  f(u),  and  its  analytical  representation: 


D(u) 


u 


u 


min 


f(u)  log^ 


f(u) 

W) 


(125) 


It  can  be  shown  that  for  x  >  0  (th's  result  may  be  obtained  from  formulas 
601.  5  and  601 . 6,  reference  19), 


(x)  2  1  -  j  . 

Combining  of  equations  125  and  126  gives 


(126) 


u  =  u 

max 

D(u)  =  ^  [f(u)] 

u  =  u 

min 


u  =  u 


min 


u  =  u 

max 


P(u) 


u  =  u 

min 


Because  of  equations  122  and  124,  the  above  inequality  becomes  D(u)  -  0. 
The  quantity  D(u)  is  equal  to  zero  only  when  for  every  value  of  u,  the  experi¬ 
mental  data,  f(u),  and  the  analytical  representation,  P(u),  are  equal.  Thus, 
the  best  representation,  P(u),  of  f(u)  minimizes  D(u): 


u  =  u 
_  max 

D(u)  =  [f(u)]  [log^  f(u)  -  log^  P(u)]  . 

u  =  u 

min 
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An  even  more  accurate  representation  of  f(u)  is  obtained  by  minimizing  the 
quantity  S(u),  which  is  a  function  of  the  squares  of  the  log  differences: 


u  =  u 


max 


S(u)  =  [f(u)]  Jlog^  f(u)  -  log^  P(u) 


u  =  u 

min 


Substitution  of  equation  123  into  equation  127  gives 


u  =  u 


max 


S(u)  =  ^  ^  [f(u)]  [log^  f{u)]  +  [f(u)]  (a  +  bu  +  cu^)  + 


u  =  u 


min 


[2f(u)  log^^  f(u)j  (a  +  bu  +  cu^) 


Minimization  of  S(u)  with  respect  to  a,  b,  and  c  provides  a  system  of 
equation  s: 

3s  as  as 

53-  -  5b  '  a^  “  °  • 

When  evaluated,  the  system  becomes 
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three 


(128) 


(129) 
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(130) 


u  =  u 


min 


The  solution  of  equations  128  through  130  provides  values  for  a,  b,  and  c. 
However,  the  solution  is  not  normalized  since  the  restraint, 


u  =  u 

max  2 

-(a  +  bu  +  cu  )  , 

e  '  -  1  , 


u  =  u 

min 

has  not  been  inserted  into  the  equation  system.  If  this  were  done,  the  resulting 
complexity  of  the  equation  system  and  its  solution  would  be  excessive.  Never¬ 
theless,  since  the  solution  is  an  approximation  of  f(u),  which  has  been  normal¬ 
ized,  it  can  be  assumed  tiiat  the  small  perturbation  in  the  value  of  a,  necessary 
to  achieve  normalization  of  P(u),  will  not  appreciably  affect  the  values  of  b  and 
c.  Therefore,  the  values  of  b  and  c  can  be  obtained  from  the  equation  system 
(equations  128  through  130)  and  the  proper  value  of  a  for  normalizing  P(u^)  be¬ 
comes 


a 


log. 


u  =  u 


E 


max 


-(bu  +  cu  ) 


u  =  u 


min 


(132) 


Consideration  of  Multicluster  Data 


In  some  cases,  the  set  of  experimental  data,  f(u),  will  exhibit  two  or  more 
distinct  regions  in  which  the  data  are  clustered.  Under  these  conditions,  the 
form  of  the  representation,  P(u),  in  equation  123  is  not  adequate.  However, 
if  the  data  set  is  divided  into  two  or  more  data  subsets  so  that  each  subset  will 
contain  only  one  cluster  of  data,  the  representation,  P(u),  in  equation  123  is 
valid.  For  example,  a  two-cluster  data  set  with  its  subsets  separated  by 
boundary  value  Ug^  is  represented  by 


P(u)  =  P^Cuj)  +  P2^^2^  ’ 


(133) 


with 


=  e 


when  u 


min 


<  u, 
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=  u 
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P,(u,)  =  0,  when  u.  >  u  , 

11  13. 


p  .  ^Z^'Z  ^  ‘^Z^Z  /  .  < 

P_{Ut)  =  e  ,  when  u  <  u_,  =  u 

Z  Z  a  Z  max 


P_,(u-)  =  0,  when  u_,  =  u  and  >  u 
Z  Z  Z  a  Z  max 

If  and  N2  represent  respectively  the  relative  portion  of  parameter  data 
points  that  lie  within  regions  and  U2  as  defined  above,  the  proper  values 
of  a^  and  a2  to  achieve  normalization  are 


=  ‘°ge  n:  ^  ‘°8, 


Cl 

2: 

u  =  u 


^  u  ^  +  c 


min 
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and 
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‘(^Z^Z  ^  ‘^Z^Z  ) 


(135) 
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APPENDIX  V 


INSTRUCTIONS  TO  SUBJECTS  FOR  TARGET -RECOGNITION  EXPERIMENTS 


HIGH -RESOLUTION  RADAR 

The  purpose  of  this  investigation  is  to  determine  whether  a  definable  relation¬ 
ship  exists  between  physical  measures  of  imagery,  such  as  brightness,  and 
human  target-recognition  performance.  The  nature  of  this  investigatic:)n  re¬ 
quires  that  many  interesting  variables,  winch  may  affei  t  performance  but  vliicii 
are  not  specifically  related  to  this  problem,  must  be  controlled  or  eliminated 
Hence,  the  task  that  you  will  be  asked  to  perform  will  bear  little,  if  any,  re¬ 
semblance  to  an  operational  situation.  The  results  of  this  scientific  study  are 
necessary  to  determine  perceptual  requirements  for  target  recognition.  Since 
the  task  in  this  investigation  is  not  similar  to  the  operational  situation,  the  re¬ 
sults  will  not  be  used  to  rate  you  in  any  way  at  your  job,  nor  will  they  be  used 
to  alter  operational  procedure. 

To  demonstrate  the  procedure  we  will  use,  1  will  initiate  the  sequence  of  op¬ 
erations  by  presenting  on  your  left-hand  display  screen  a  picture  of  the  target 
you  will  be  asked  to  find.  This  picture  will  be  the  same  size  and  in  the  san  e 
orientation  as  the  target  in  the  radar  imagery.  You  will  study  this  cue  as  long 
as  you  feel  necessary  and  inform  me  when  you  are  ready.  1  will  then  present 
the  high  -  re  solution  radar  imagery  on  the  center  screen.  Some  targets  will  be 
very  easy  to  find,  others  will  be  extremely  difficult.  Find  the  target  that  you 
have  studied  as  soon  as  ycm  can.  We  are  interested  in  both  the  speed  and  ac¬ 
curacy  with  which  you  can  locate  these  targets.  When  you  locate  the  target, 
press  the  top  switch  on  your  left.  This  action  will  cause  a  magnified  version 
of  the  imagery  right  around  the  cross  hairs  to  appear  on  the  right-hand  screen. 
Then  move  the  cross  hairs  over  the  target.  When  you  have  laid  the  cross  hairs, 
press  the  second  switch  on  your  left.  If  at  some  point  you  realize  tliat  you  have 
erred  in  your  identification,  inform  me,  search  again  for  the  target,  and  lav 
the  cross  hairs  over  it.  We  will  now  have  several  practice  trials  to  fam.iliarize 
you  with  the  equipment  and  procedure. 

INFRARED  AND  OPTICAL  IMAGERY 

The  purpose  of  our  investigation  is  to  determine  if  some  relationship  can  be 
found  between  physical  measures  of  imagery,  such  as  brightness  or  variations 
in  brightness  between  a  target  and  its  surrounding  area,  and  human  perform¬ 
ance  in  a  target  recognition  task.  That  is,  is  a  target  that  is  very  bright  in  a 
background  of  moderiite  brightness  easier  or  more  difficult  to  locate  than  a 
similar  target  that  is  very  dark?  Does  the  difference  between  a  target  and  its 
background  necessary  for  the  former  to  be  discernible  vary  with  the  size  of 
the  target?  We  are  ti  ying  to  determine  relationships  of  this  nature  using  vari¬ 
ous  types  of  imagery.  In  this  particular  phase  of  the  investigation,  we  are 
using  optical  and  infrared  imagery. 

To  demonstrate  the  procedure,  I  will  first  show  you  a  picture  of  the  target 
you  are  to  find.  This  is  an  exact  copy  of  the  target,  with  one  exception.  The 
transparent  tape  over  the  target  makes  the  target  key  appear  slightly  darker 
than  it  will  appear  on  the  imagery.  The  target  key  will  also  be  in  the  same 
orientation  as  tlie  target  in  the  imagery.  1  will  tell  you  whether  you  will  be 
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luokiii^  lit  infrared  nr  optical  photography  and  what  c  ategory  the  target  is  in; 
tliat  IS,  wlu-ther  yon  will  be  looking  for  a  shopping  center,  industrial  complc-x, 
etc  .  Study  the  target  key  as  long  as  you  wish.  When  you  are  ready,  I  will 
turn  over  the  photograph  in  which  you  are  to  locate  the  targc’t.  Find  the  tar¬ 
get  as  quickly  as  you  can,  but  be  sure  to  be  accurate.  I  will  measure  the  time 
you  take  to  locate  the  target. 

I  have  si.\  demonstration  photographs  to  better  acquaint  you  with  the  procedure 
and  task.  Ask  any  questions  you  like  bc'fore  we  go  on  with  the  test  photographs. 
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